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Introduction: 

Our  goal  is  to  develop  a  diagnostic  device  for  breast  cancer  based  on  fluorescence 
and  diffuse  reflectance  spectroscopy.  Our  specific  objectives  are  to  characterize  the 
fluorescence  properties  of  normal,  benign,  and  malignant  breast  tissue  in  patients 
undergoing  breast  cancer  and  breast  reduction  surgeries.  Additionally,  we  will 
characterize  the  intrinsic  fluorescence  properties  of  tissue  constituents  in  order  to  identify 
the  biological  basis  for  spectral  differences  in  tissue.  An  initial,  ex  vivo,  study  will  be 
completed,  and  then  a  fast,  portable  device  will  be  employed  to  measure  the  fluorescence 
and  diffuse  reflectance  properties  of  tissue,  in  vivo.  Results  will  be  correlated  with  the 
gold  standard,  histopathology,  to  develop  statistical  methods  to  determine  the  diagnostic 
potential  of  this  technology. 

Body: 

The  first  major  focus  of  work  over  the  past  12  months  has  been  in  identifying  the 
biological  basis  for  the  effectiveness  of  fluorescence  spectroscopy  in  diagnosing  breast 
cancer.  To  this  end,  the  fluorescence  properties  of  one  normal  and  four  malignant  cell 
lines  were  characterized  using  fluorescence  spectroscopy.  Specifically,  the  fluorescence 
properties  of  three  intrinsic  fluorophores:  tryptophan,  nicotinamide  adenine  dinucleotide 
(NAD(P)H),  and  flavin  adenine  dinucleotide  (FAD),  were  characterized  using 
fluorescence  spectroscopy.  Additionally,  two-photon  excitation  microscopy  was 
employed  to  image  NAD(P)H  fluorescence  from  the  normal  and  one  of  the  malignant  cell 
lines.  It  was  found  that  the  fluorescence  intensity  of  tryptophan  and  NAD(P)H  were 
significantly  reduced  in  the  malignant  cell  lines  (p  <  0.05).  A  manuscript  detailing  these 
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results  and  their  implications  has  been  accepted  into  the  peer-reviewed  journal, 
Photochemistry  and  Photobiology,  and  is  included  in  Appendix  A  (p.6).  This  work 
addresses  part  (a)  of  Task  1  in  the  Statement  of  Work. 

Parts  (b-c)  of  Task  1  in  the  statement  of  work  involve  the  measurement  of 
fluorescence  and  diffuse  reflectance  spectra  of  human  breast  tissues,  ex  vivo. 
Unfortunately,  some  difficulty  was  encountered  in  gaining  approval  for  conducting 
human  subjects  research  from  the  DOD  Institutional  Review  Board.  Approval  has 
recently  been  obtained,  however,  and  we  plan  to  proceed  with  collection  of  this  data  in 
the  coming  year.  This  has  delayed  the  start  of  our  study,  but  should  not  affect  our  overall 
goals. 

A  preliminary  study  involving  another  protocol  not  covered  under  this  grant  was 
completed  however,  and  serves  to  answer  some  of  the  questions  posed  by  Task  1  of  the 
Statement  of  Work.  A  comparison  of  classification  accuracy  using  empirically -based 
statistical  methods  was  carried  out.  Fluorescence  and  diffuse  reflectance  properties  of 
breast  tissue  was  characterized  for  a  total  of  56  samples  from  32  patients,  immediately 
after  surgery.  To  briefly  summarize,  it  was  found  that  fluorescence-based  detection  of 
breast  cancer  was  highly  successful,  with  a  cross-validated  classification  accuracy  of 
84%.  It  was  found  that  diffuse  reflectance-based  detection  of  breast  cancer  was  less 
effective,  with  a  cross  validated  classification  accuracy  of  61%.  A  manuscript  detailing 
the  results  of  this  study  was  accepted  for  publication  into  the  peer-reviewed  journal,  IEEE 
Transactions  on  Biomedical  Engineering  and  is  included  in  Appendix  B  (p.  14). 

One  limitation  of  this  study  is  that  it  employed  a  fixed  probe  geometry,  with  a 
central  bundle  of  optical  fibers  for  collection  of  light,  and  an  outer  ring  of  illumination 
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fibers.  Furthermore,  little  research  has  been  done  to  identify  how  probe  geometry  affects 
the  diagnostic  efficacy  of  these  methods.  The  study  to  be  completed  under  this  grant’s 
protocol  will  be  carried  out  using  a  probe  designed  with  multiple  separation  distances 
between  the  illumination  and  collection  fibers.  A  schematic  of  this  probe  design  is 
shown  in  Fig.  1.  Because  the  probe  geometry  affects  the  probing  volume  within  the 
tissue  [1],  having  multiple  probing  geometries  will  allow  us  to  evaluate  in  what  way 
probe  geometry  affects  the  information  content  obtained  with  these  measurements,  and 
will  also  help  determine  if  there  is  an  optimal  probe  geometry  for  use  in  future  studies. 


Fig.  1 :  Probe  geometry  to  be  used  in  the  current  human  subjects  protocol.  The  probe  consists  of 
a  central  illumination  region  (black)  surrounded  by  three  concentric  rings  of  collection  fibers. 

Part  (d)  of  Task  1  pertains  to  developing  a  method  of  analysis  for  interpreting 
spectral  data  taken  from  the  tissue  specimens.  To  this  end  a  novel  statistical  method  of 
analysis,  appropriate  for  small  sizes,  has  been  developed  and  tested  on  spectral  data 
obtained  from  the  preliminary  ex  vivo  breast  tissue  study.  This  method  is  described  in 
detail  in  the  manuscript  submitted  in  Appendix  B  (p.  14),  and  will  be  applied  to  future 
data  collected  in  this  study. 

Key  Research  Accomplishments 

•  Identification  of  intrinsic  differences  in  fluorescence  properties  of  human  breast 
cell  cultures  with  malignancy 
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•  Identification  of  spectral  differences  present  in  human  breast  tissue,  ex  vivo. 

•  Development  of  a  statistical  method  of  analysis  for  evaluating  the  classification 
accuracy  of  optical  spectroscopy  in  diagnosing  breast  cancer 

Reportable  Outcomes: 

Degrees  Obtained: 

Master’s  of  Science  -  Biomedical  Engineering,  Dec.  2002 
Journal  Articles: 

•  Palmer,  G.M.,  C.  Zhu,  T.M.  Breslin,  F.  Xu,  K.W.  Gilchrist,  and  N.  Ramanujam. 
“Comparison  of  multi-excitation  fluorescence  and  diffuse  reflectance 
spectroscopy  for  the  diagnosis  of  breast  cancer,  ex  vivo.”  Accepted,  IEEE  Trans 
BME. 

•  Breslin,  T.M.,  F.  Xu,  G.M.  Palmer,  C.  Zhu,  K.W.  Gilchrist,  and  N.  Ramanujam. 
“Auto  fluorescence  and  diffuse  reflectance  properties  of  malignant  and  benign 
breast  tissues.”  Accepted,  Ann.  Surg.  One. 

•  Palmer,  G.M.,  P.J.  Keely,  T.M.  Breslin,  and  N.  Ramanujam.  “Autofluorescence 
spectroscopy  of  normal  and  malignant  human  breast  cell  lines.”  Accepted, 
Photochem  Photobiol. 

Conference  Presentations: 

•  Palmer,  G.M.,  C.  Zhu,  T.M.  Breslin,  P.J.  Keely,  F.  Xu,  K.W.  Gilchrist,  J.M. 
Squirrell,  J.G.  White,  and  N.  Ramanujam.  “Use  of  auto  fluorescence  and  diffuse 
reflectance  spectroscopy  for  the  diagnosis  of  breast  cancer,  ex  vivo,  and  its 
biological  basis.”  Poster  Presentation,  Third  Inter-Institute  Meeting  on 
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Diagnostic  Optical  Imaging  and  Spectroscopy:  The  Clinical  Adventure.  National 
Institutes  of  Health,  Bethesda,  MD,  Sept.  2002. 

•  Palmer,  G.M.,  P.J.  Keely,  T.M.  Breslin,  J.M.  Squirrell,  J.G.  White,  and  N. 
Ramanujam.  “Understanding  Sources  of  Contrast  in  Fluorescence  Based 
Diagnosis  of  Breast  Cancer.”  Poster  Presentation,  Engineering  Conferences 
International:  Advances  in  Optics.  Banff,  Canada.  Aug.  2003. 

Conclusions 

Significant  differences  in  the  fluorescence  and  diffuse  reflectance  properties  of 
normal  and  malignant  breast  tissues  have  been  identified,  and  a  method  for  classification 
of  these  tissues  has  been  developed.  A  diagnostic  device  based  on  fluorescence  and 
diffuse  reflectance  spectroscopy  has  the  advantage  of  being  fast,  quantitative,  and 
minimally  invasive,  and  could  improve  breast  cancer  care  both  in  assessing  surgical 
margins  for  residual  cancer  in  vivo  during  surgery,  as  well  as  during  core  needle  biopsy  to 
accurately  identify  the  appropriate  region  to  biopsy.  Furthermore,  the  biological  basis  for 
observed  differences  in  fluorescence  properties  of  the  tissue  has  been  examined,  and 
differences  in  the  intrinsic  fluorescence  properties  of  normal  and  malignant  human  breast 
cell  lines  have  been  identified.  These  differences  in  intrinsic  cellular  fluorescence  could 
explain,  in  part,  the  biological  basis  for  observed  spectral  differences  in  tissue. 
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ABSTRACT 

The  fluorescence  of  tryptophan,  reduced  nicotinamide  adenine 
dinucleotide  (phosphate)  (NAD(P)H)  and  flavin  adenine  di¬ 
nucleotide  (FAD)  were  characterized  in  normal  human  breast 
cells  as  well  as  in  malignant  human  breast  cells  of  similar  and 
dissimilar  genetic  origins.  Fluorescence  measurements  of  each 
cell  line  were  made  over  a  wide  range  of  cell  concentrations, 
and  the  fluorescence  per  cell  was  determined  from  the  slope 
in  the  linear  range  of  the  fluorescence  intensity  cell  con¬ 
centration  plot.  All  of  the  malignant  cells  showed  a  statistically 
signiflcant  decrease  in  the  tryptophan  fluorescence  per  cell 
relative  to  that  of  the  normal  cells.  No  statistically  signiflcant 
differences  were  observed  in  the  NAD(P)H  or  FAD  fluores¬ 
cence  per  cell  between  the  normal  and  any  of  the  malignant 
cell  types.  NAD(P)H  fluorescence  was  also  imaged  from 
monolayers  of  the  normal  and  malignant  cells  (of  similar 
genetic  origin)  using  two-photon  fluorescence  microscopy.  A 
statistically  signiflcant  decrease  in  the  NAD(P)H  fluorescence 
with  malignancy  was  observed,  suggesting  that  fluorescence 
imaging  of  single  cells  or  the  cell  monolayer  preparation  may 
provide  more  contrast  than  volume-averaged  fluorescence 
measurements  of  cells  in  suspension.  In  conclusion,  the  dif¬ 
ferences  in  normal  and  malignant  human  breast  tissue  fluo¬ 
rescence  spectra  may  be  attributed  in  part  to  differences  in  the 
intrinsic  cellular  fluorescence  of  normal  and  malignant  breast 
epithelial  cells. 


INTRODUCTION 

Breast  cancer  screening  relies  on  mammography,  which  is  effective 
in  identifying  suspicious  lesions.  However,  a  significant  number  of 
false  positives  are  produced  (1).  As  a  result,  follow-up  diagnostic 
procedures  such  as  percutaneous  core  needle  biopsy  with  image 
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guidance  are  performed  to  further  evaluate  these  abnormalities  (1). 
Core  needle  biopsy  has  a  limited  sampling  accuracy  because  only 
a  few  small  pieces  of  tissue  are  extracted  from  random  locations  in 
the  suspicious  mass.  In  some  cases,  the  suspicious  mass  may  be 
missed  altogether. 

Optical  sensors  based  on  ultraviolet-visible  (UV-VIS)  spec¬ 
troscopy  have  the  potential  to  improve  the  sampling  accuracy  of 
core  needle  biopsy.  Optical  techniques  have  the  capability  to 
survey  multiple  tissue  sites  without  the  need  for  biopsy  and  provide 
real-time  diagnosis.  UV-VIS  fluorescence  spectroscopy  is  a  tech¬ 
nique  that  has  been  shown  to  have  consistently  high  sensitivities 
and  specificities  for  detecting  human  precancers  and  early  cancers 
in  vivo  in  a  number  of  organ  sites,  including  the  gastrointestinal 
tract,  bronchus  and  cervix  (2).  Incorporating  UV-VIS  fluorescence 
spectroscopy  as  an  adjunct  diagnostic  modality  to  image-guided 
core  needle  biopsy  could  potentially  reduce  the  false-negative  rate 
of  current  breast  biopsy  procedures.  This  in  turn  could  lead  to 
fewer  follow-up  procedures  and  fewer  repeat  biopsies  in  patients 
suspected  to  have  breast  cancer.  A  second  potential  application  of 
UV-VIS  fluorescence  spectroscopy  is  the  evaluation  of  tumor 
margins  during  breast  cancer  surgery.  A  real-time,  accurate  means 
of  assessing  tumor  margins  could  reduce  intraoperative  time  as 
well  as  provide  greater  assurance  that  the  entire  lesion  has  been 
removed. 

Breast  tissue,  which  is  primarily  composed  of  epithelial  cells,  an 
extracellular  matrix  and  fat,  contains  a  number  of  endogenous 
fluorophores.  Tryptophan,  reduced  nicotinamide  adenine  dinucle¬ 
otide  (phosphate)  (NAD(P)H)  and  flavin  adenine  dinucleotide 
(FAD)  are  present  in  the  epithelial  cells  (3-6).  The  excitation- 
emission  maxima  corresponding  to  the  peaks  of  these  fluorophores 
are  280-340  nm  for  tryptophan,  340-460  nm  for  NAD(P)H  and 
460-520  nm  for  FAD  (3-6).  Tryptophan  accounts  for  the  majority 
of  protein  fluorescence,  and  its  emission  is  sensitive  to  the  polarity 
of  the  environment  (7).  NAD(P)H  and  FAD  are  involved  in  the 
oxidation  of  fuel  molecules  and  can  be  used  to  probe  changes  in 
cellular  metabolism  (8).  The  primary  fluorophore  in  the  breast 
tissue  extracellular  matrix  is  type  I  collagen  (9),  which  has  a  peak 
near  330-400  nm  that  is  attributed  to  cross-links  (4).  Other  major 
factors  that  influence  breast  tissue  fluorescence  spectra  in  the  UV- 
VIS  spectrum  are  endogenous  absorbers  and  scatterers  in  the  tissue. 
Absorption  in  tissue  is  dominated  by  oxygenated  and  deoxygen- 
ated  hemoglobin  (10).  Scattering  occurs  because  of  local  changes 
in  the  index  of  refraction  across  mitochondrial  and  other  membrane 
interfaces  (10).  Wavelength-dependent  absorption  and  scattering 
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affect  both  the  spectral  intensity  and  the  line  shape  of  the  fluo¬ 
rescence  measured  from  the  tissue. 

Several  groups  have  demonstrated  that  there  are  significant  dif¬ 
ferences  in  the  fluorescence  excitation  and  emission  spectra  of 
normal,  benign  and  malignant  breast  tissues  ex  vivo  (1 1-17).  How¬ 
ever,  it  is  difficult  to  directly  interpret  the  sources  of  endogenous 
fluorescence  contrast  between  the  malignant  and  nonmalignant 
breast  tissues  from  the  measured  tissue  fluorescence  spectra.  This  is 
due  to  the  complex  interplay  of  fluorescence,  absorption  and  scat¬ 
tering.  Thus,  systematic  investigations  of  breast  tissue  constituent 
fluorescence  in  progressively  complex  biological  systems  that 
range  from  cells  to  tissue  cultures  to  living  animal  models  of  breast 
cancer  may  provide  an  effective  approach  toward  understanding 
the  differences  in  the  endogenous  fluorescence  of  malignant  and 
nonmalignant  breast  tissues. 

Several  groups  have  characterized  the  fluorescence  of  normal 
and  malignant  cells  for  a  number  of  cell  types  (5,18-22).  Thus  far, 
only  one  group  has  studied  the  fluorescence  properties  of  normal 
and  malignant  human  breast  cells  (22).  Glassman  et  al.  (22) 
measured  fluorescence  from  the  normal  HTB125  and  the  malignant 
HTB22  and  HTB126  human  breast  cells  (American  Type  Culture 
Collection,  Rockville,  MD).  The  HTB125  and  HTB126  were  de¬ 
rived  from  the  same  patient.  Fluorescence  spectra  were  recorded 
at  an  excitation  wavelength  of  353  nm  from  cells  that  were 
centrifuged  and  then  packed  into  a  cuvette.  The  results  showed  that 
the  NAD(P)H  fluorescence  of  the  malignant  cells  is  a  factor  of  2 
higher  than  that  of  normal  cells  at  this  excitation  wavelength.  The 
observed  increase  in  NAD(P)H  fluorescence  was  attributed  to  the 
increased  metabolic  rate  of  the  malignant  cells. 

These  results  set  an  important  precedent  for  future  work. 
However,  there  are  a  number  of  potential  refinements  that  can  be 
incorporated  into  future  studies  that  will  improve  the  quality  and 
information  content  of  cellular  fluorescence  measurements.  This 
includes  (1)  the  comparison  of  normal  and  malignant  cell  phe¬ 
notypes  of  genetically  similar  origins;  (2)  characterization  of  all  the 
endogenous  fluorophores  present  in  these  cells;  (3)  measurement  of 
cellular  fluorescence  at  concentrations  at  which  the  effects  of  tur¬ 
bidity  are  minimal;  and  (4)  viability  testing  before  and  after  the 
experimental  measurements. 

The  goal  of  this  work  was  to  systematically  characterize  the 
differences  in  the  fluorescence  properties  of  normal  and  malignant 
human  breast  cell  lines.  This  study  builds  upon  previous  work  and 
attempts  to  minimize  potential  sources  of  error  that  can  influence 
the  analysis  and  interpretation  of  fluorescence  measurements  from 
cells.  Measurements  were  made  from  suspensions  of  breast  cells  to 
allow  for  the  characterization  of  fluorescence  from  a  relatively 
simple  biological  system  of  the  normal  breast  and  breast  cancer 
cells.  In  this  system,  there  was  no  direct  cell  to  cell  or  cell  to 
extracellular  matrix  interaction.  These  effects,  which  are  important 
for  cell  signaling,  are  evaluated  to  a  limited  extent  and  will  be 
evaluated  more  thoroughly  in  an  independent  study  that  includes 
cell  monolayers  and  tissue  cultures  (see  Discussion).  These  step¬ 
wise  investigations  will  enable  a  systematic  comparison  of  the 
effects  of  different  microenvironments  on  the  fluorescence  prop¬ 
erties  of  normal  and  malignant  human  breast  cells. 

Three  cell  lines  were  used  in  these  studies.  These  included 
the  well-differentiated  T47D  cells  and  the  poorly  differentiated 
MDA231  cells,  which  are  both  spontaneously  occurring  human 
breast  adenocarcinomas,  and  the  MCFIO  cell  line,  which  is  a 
putatively  normal  human  breast  cell  line.  The  T47D  cell  line  was 
derived  from  a  54-year-old  patient  with  infiltrating  ductal 


carcinoma  (23).  The  MDA231  cell  line  was  derived  from  a  51- 
year-old  woman  with  cancer  that  had  developed  a  plural  effusion 
(24).  Finally,  the  MCFIO  cell  line  was  derived  from  benign 
mastectomy  tissue  obtained  from  a  36-year-old  woman  (25).  These 
cells  became  spontaneously  immortalized  in  culture  into  a  non- 
tumorigenic,  putatively  normal  cell  line  (25).  To  compare  normal 
and  malignant  phenotypes  of  similar  genetic  origins,  MCFIO  cells 
were  transformed  by  stable  expression  of  the  R-Ras  oncogene 
(malignant  phenotype)  or  transfected  with  only  the  pZIP  vector 
(normal  phenotype).  Although  the  immortalized  MCFIO  cells  may 
not  be  “truly  normal,”  transforming  these  cells  with  the  R-Ras 
oncogene  enables  a  direct  comparison  between  transformed  cells 
(with  a  known  oncogene)  and  their  nontransformed  counterpart.  In 
addition,  the  T47D  cell  line  was  stably  transfected  with  the  R-Ras 
oncogene  to  serve  as  a  positive  control  because  a  normal  to  malig¬ 
nant  transformation  is  not  expected  in  these  cells. 

Ruorescence  spectra  of  suspensions  from  each  cell  type  (from 
a  total  of  five  cell  types)  were  measured  over  the  entire  UV-VIS 
range  to  characterize  the  fluorescence  of  all  three  endogenous 
fluorophores,  tryptophan,  NAD(P)H  and  FAD,  present  in  the  cells. 
For  each  cell  type,  the  fluorescence  measurements  were  carried  out 
over  a  range  of  cell  concentrations,  and  the  fluorescence  per  cell 
was  determined  from  the  slope  in  the  linear  range  of  the  fluo¬ 
rescence  intensity  vs  cell  concentration  plot.  Cell  viability  was 
determined  before  and  after  each  experiment.  In  addition,  the  cell 
sizes  were  quantitatively  determined  for  each  cell  type. 

All  of  the  malignant  cells  showed  a  statistically  significant 
decrease  in  the  tryptophan  fluorescence  per  cell  relative  to  that  of 
the  normal  cells.  No  statistically  significant  differences  were 
observed  in  the  NAD(P)H  or  FAD  fluorescence  per  cell  between 
the  normal  and  any  of  the  malignant  cell  types.  The  tryptophan 
fluorescence  per  ceU  was  found  to  be  very  similar  for  the  T47D  and 
T47D  R-Ras  (positive  control)  ceU  lines,  confirming  that  the 
differences  in  the  tryptophan  fluorescence  of  the  nontransformed 
MCFIO  pZIP  cells  and  the  transformed  MCFIO  R-Ras  cells  are  not 
due  to  the  expression  of  the  R-Ras  oncogene  but  rather  due  to  the 
transformation  of  the  MCFIO  cells  from  a  normal  to  a  malignant 
phenotype.  In  addition,  cell  viability  (greater  than  80%)  and  cell 
sizes  (150-200  pm^)  were  found  to  be  similar  for  all  cell  lines 
evaluated  in  this  study. 

NAD(P)H  fluorescence  was  also  imaged  from  monolayers  of  the 
normal  and  malignant  cells  (of  similar  genetic  origin)  using  two- 
photon  fluorescence  microscopy.  A  statistically  significant  decrease 
in  the  NAD(P)H  fluorescence  with  malignancy  was  observed, 
suggesting  that  fluorescence  imaging  of  single  cells  or  the  cell 
monolayer  preparation  may  provide  more  contrast  than  volume- 
averaged  fluorescence  measurements  of  cells  in  suspension. 

MATERIALS  AND  METHODS 

Preparation  of  cell  suspensions.  Cell  lines  were  obtained  from  the 
American  Type  Culture  Collection  and  were  maintained  free  of  myco¬ 
plasma  or  other  contamination  by  adherent  culture  according  to  their  indi¬ 
vidual  established  protocols.  Measurements  were  made  on  cells  having 
passage  numbers  of  50  or  less.  The  cells  were  not  synchronized.  The  T47D 
cells  were  grown  in  Roswell  Park  Memorial  Institute  1640  medium 
supplemented  with  10%  fetal  bovine  serum  and  8  pg/mL  insulin.  The 
MDA231  cells  were  grown  in  Dulbecco  modified  Eagle  media  (DMEM)- 
F12  supplemented  with  10%  horse  serum.  The  MCFIO  cells  were  grown  in 
DMEM-F12  supplemented  with  5%  horse  serum,  10  pg/mL  insulin,  20  ng/ 
mL  epidermal  growth  factor  and  0.5  pg/mL  hydrocortisone.  Cells  were 
incubated  at  3TC  with  5%  CO2  for  the  T47D  cells  and  10%  CO2  for  the 
MDA231  and  MCFIO  cells.  Transfection  of  the  MCFIO  and  T47D  cells 
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with  the  R-Ras  oncogene  was  performed  according  to  previously  described 
protocols  (26). 

Upon  reaching  confluence  (generally  3-4  days  after  plating),  the  cells 
were  detached  from  their  plates  by  trypsinization  and  triply  washed  in 
phosphate-buffered  saline  (PBS).  This  was  carried  out  using  10  mL  of  PBS 
per  wash.  After  the  final  wash,  the  cells  were  resuspended  in  a  final  volume 
of  3  mL  PBS.  Cell  concentration  was  determined  by  counting  the  number  of 
cells  per  milliliter  using  a  hemacytometer.  This  measurement  was  repeated 
four  times  (by  making  an  independent  count  in  each  of  the  four  grids  in  the 
hemacytometer)  to  obtain  a  standard  deviation  of  the  concentration  of  cells. 
The  solution  containing  the  cells  was  then  diluted  to  attain  the  experimental 
concentration.  Cell  suspensions  at  approximately  eight  concentrations  (5  X 
10^  cells/mL  to  3  X  10^  cells/mL),  each  taken  from  an  independent  cell 
culture,  were  characterized  for  each  cell  line. 

The  viability  was  visually  evaluated  by  checking  membrane  integrity 
using  a  phase-contrast  stereomicroscope.  This  process  was  repeated  four 
times  (again  using  each  of  the  four  hemacytometer  grids)  to  obtain 
a  standard  deviation  for  cell  viability.  The  accuracy  of  determining  cell 
viability  by  sight  with  a  phase-contrast  microscope  was  compared  with 
a  trypan  blue  exclusion  assay  (gold  standard).  The  viability  assays  using 
phase-contrast  microscopy  and  the  trypan  blue  assay  were  carried  out  12 
times  (four  repetitions  at  each  of  three  different  concentrations)  from  a  T47D 
cell  culture  over  a  range  of  concentrations  within  the  experimental 
concentration  range  (5X10^  cells/mL  to  3  X  10^  cells/mL).  The  average 
viability  determined  by  phase-contrast  microscopy  and  the  trypan  blue 
assay  was  87.5  ±  5.9%  and  88.1  ±  7.6%,  respectively.  A  Wilcoxon  rank- 
sum  test  (27)  was  used  to  confirm  that  the  differences  in  cell  viability 
obtained  with  the  two  methods  are  statistically  insignificant  (F  >  0.5). 

Fluorescence  spectroscopy  of  cell  suspensions.  The  cell  suspension  was 
placed  in  a  1  cm  path  length  quartz  cuvette.  A  convention^  laboratory 
fluorometer  (Fluorolog-3,  JY  Horiba  Inc.,  Edison,  NJ)  was  used  to  char¬ 
acterize  the  fluorescence  properties  of  each  cell  suspension.  The  instrument 
consists  of  a  450  W  xenon  arc  lamp,  double  excitation  and  emission 
scanning  monochromators,  a  sample  compartment  with  a  cuvette  holder  and 
magnetic  stirrer  and  a  photomultiplier  tube  (PMT).  Adjustable  parameters 
include  the  excitation  and  emission  wavelength  ranges  and  increment,  the 
excitation  and  emission  band-passes,  the  high  voltage  of  the  PMT,  the 
integration  time,  the  stirring  speed  of  the  magnetic  stirrer  and  the  collection 
geometry  (right  angle  or  front  face). 

Fluorescence  excitation-emission  matrices  (EEM)  were  recorded  by 
measuring  the  fluorescence  emission  spectra  over  a  range  of  excitation 
wavelengths.  A  fluorescence  EEM  involves  measurement  of  the  fluores¬ 
cence  intensity  as  a  function  of  both  excitation  and  emission  wavelengths. 
The  excitation  wavelength  was  varied  from  260  to  540  nm  in  20  nm 
increments.  At  each  excitation  wavelength,  the  fluorescence  emission  spec¬ 
tra  were  recorded  from  275  to  700  nm  in  5  nm  increments.  These 
wavelength  ranges  enabled  characterization  of  all  endogenous  fluorophores, 
tryptophan,  NAD(P)H  and  FAD,  present  in  cells  in  the  UV-VIS  spectrum. 
The  excitation  and  emission  band-passes  were  2  and  2.5  nm,  respectively. 
The  PMT  high  voltage  was  set  at  950  V,  and  an  integration  time  of  1  s  per 
wavelength  was  used  in  all  experiments  to  achieve  an  adequate  signal-to- 
noise  ratio  while  minimizing  the  measurement  time.  This  resulted  in  a  total 
measurement  time  of  30  min  for  each  cell  suspension.  A  magnetic  stirrer 
was  used  to  keep  the  cells  in  suspension  during  the  course  of  the  30  min 
measurement  time.  The  stirring  speed  was  determined  by  finding  the 
minimum  speed  at  which  fluorescence  intensity  remained  constant  over 
time  to  minimize  any  potential  physical  damage  to  the  cells.  Right-angle 
collection  was  used  to  measure  the  fluorescence  EEM  from  each  cell 
suspension  in  order  to  minimize  the  detection  of  the  backscattered  exci¬ 
tation  light.  After  each  experiment,  a  repeat  fluorescence  spectrum  was 
measured  at  260  nm  excitation  and  compared  with  that  measured  at  the 
beginning  of  the  experiment  to  establish  that  the  measurement  protocol  used 
did  not  induce  photobleaching.  The  peak  intensity  at  an  excitation-emission 
wavelength  of  260-340  nm  of  the  repeat  fluorescence  spectrum  was  found 
to  be  within  10%  of  that  taken  at  the  beginning  of  the  experiment.  In 
addition,  the  fluorescence  intensity  at  the  excitation-emission  maxima  in 
the  EEM  corresponding  to  tryptophan,  NAD(P)H  and  FAD  was  measured 
once  per  minute  (Is  integration  time)  for  a  total  of  5  min.  The  average  and 
standard  deviation  of  the  five  data  points  provided  the  reproducibility  of  the 
fluorescence  measurements  obtained  from  each  cell  suspension. 

Fluorescence  measurements  were  obtained  from  the  cell  suspension  as 
well  as  the  PBS  media;  the  latter  was  subtracted  as  background  fluo¬ 
rescence.  The  fluorescence  intensity  at  each  wavelength  of  the  emission 
spectrum  (for  a  particular  excitation  wavelength)  was  divided  by  that 


measured  with  a  reference  photodiode  in  the  excitation  monochromator, 
this  correction  accounts  for  the  wavelength-dependent  variations  in  the 
lamp  intensity  and  the  wavelength-dependent  efficiency  of  the  excitation 
monochromator.  A  factory-generated  correction  factor  was  also  used  to 
correct  for  the  wavelength-dependent  efficiency  of  the  optical  components 
in  the  emission  monochromator  and  the  PMT.  Finally,  the  fluorescence 
intensity  at  each  excitation-emission  wavelength  pair  measured  from  the 
cells  was  divided  by  the  peak  fluorescence  intensity  (excitation-emission: 
460-575  nm)  of  a  Rhodamine  standard  dissolved  in  water  (2.14  pM)  to 
correct  for  day-to-day  variations  in  the  throughput  of  the  instrument. 

Data  analysis  to  determine  the  fluorescence  intensity  per  cell.  A  Unear 
relationship  is  expected  between  fluorescence  intensity  and  cell  concentra¬ 
tion  when  turbidity  is  negligible,  and  the  slope  of  the  linear  relationship 
represents  the  fluorescence  intensity  per  cell.  Thus,  the  fluorescence  in¬ 
tensity  at  each  excitation-emission  maximum  was  plotted  as  a  function  of 
the  concentration  of  cells  in  suspension.  The  excitation-emission  maxima 
of  tryptophan,  NAD(P)H  and  FAD  are  280-340  nm,  340-460  nm  and  460- 
520  nm,  respectively.  Although,  tryptophan  was  observed  to  maximally 
fluoresce  at  280-340  nm,  fluorescence  measurements  taken  at  several  of  the 
higher  cell  concentrations  saturated  the  PMT.  Thus,  the  tryptophan  fluo¬ 
rescence  measured  from  the  cells  was  analyzed  at  an  excitation-emission 
wavelength  of  300-340  nm. 

A  Unear  range  was  quantitatively  determined  by  doing  a  series  of  linear 
least  squares  fits  (27),  in  which  data  points  at  the  upper  range  of  con¬ 
centrations  were  sequentially  excluded  until  the  value  was  maximized. 
Because  the  intercept  was  expected  to  pass  through  the  origin,  it  was  fixed 
at  zero  for  each  fit.  However,  it  should  be  noted  that  similar  results  were 
obtained  whether  or  not  the  intercept  of  the  fit  was  fixed  at  zero.  The  slope 
and  the  standard  errors  of  the  slope  were  calculated  using  the  R  software 
package  (27).  The  slope  reflects  the  fluorescence  intensity  per  cell  and 
enables  a  direct  comparison  of  the  intrinsic  fluorescence  between  cell  lines. 
Because  each  data  point  in  the  plot  was  obtained  from  an  independent 
culture,  the  standard  error  of  the  slope  reflects  the  interculture  variation  in 
the  fluorescence  measurement. 

Determination  of  cell  sizes.  Variations  in  ceU  size  may  affect  the  amount 
of  fluorescence  measured  per  cell.  Therefore,  the  average  size  of  the  cells 
was  determined  for  each  cell  line.  The  cells  were  imaged  using  a  phase- 
contrast  microscope  coupled  to  a  CCD  camera.  These  images  were  impor¬ 
ted  into  Image!  (NIH,  http://rsb.info.nih.gov/ij).  Cell  areas  were  determined 
using  the  “adjust  threshold”  and  “analyze  particles”  functions.  These  func¬ 
tions  identify  areas  contained  within  the  cells  and  thus  provide  an  area  of 
distribution.  This  was  done  for  approximately  20  cells  within  an  inde¬ 
pendent  culture  for  each  cell  type.  Tfte  microscope  was  calibrated  such  that 
cell  sizes  could  be  reported  in  absolute  units. 

Two-photon  fluorescence  microscopy.  The  NAD(P)H  fluorescence  in 
monolayers  of  MCFIO  pZIP  and  MCFIO  R-Ras  cells  was  imaged  using 
two-photon  fluorescence  microscopy.  These  experiments  were  carried  out 
to  compare  the  fluorescence  measurements  obtained  from  cell  monolayers 
with  those  obtained  from  cells  in  suspension.  MCFIO  pZIP  cells  and 
MCFIO  R-Ras  cells  were  plated  separately  on  lysine-coated  cover  glass 
bottom  petri  dishes,  immersed  in  their  standard  growth  media  and  allowed 
to  grow  overnight  under  the  same  conditions  as  in  the  previous  experiment 
(see  Preparation  of  Cell  Suspensions).  Three  pairs  of  MCFIO  pZIP  and 
MCFIO  R-Ras  cultures  were  imaged  using  two-photon  fluorescence  mi¬ 
croscopy.  Both  cell  lines  within  each  pair  were  imaged  under  identical 
experimental  conditions. 

The  two-photon  excitation  source  was  a  Ti:  Sapphire  laser  tuned  to  either 
775  nm  or  789  nm  (this  difference  was  due  to  day-to-day  variations  in  the 
stability  of  the  laser).  The  excitation  light  was  delivered  to  and  the  emitted 
light  was  collected  from  the  sample  through  a  conventional  inverted  Nikon 
eclipse  microscope  (with  a  40x  oil  immersion  objective  with  a  numerical 
aperture  of  1 .3)  coupled  to  a  Biorad  confocal  scan  head.  A  PMT  was  used 
as  the  detector.  The  absolute  light  intensity  reaching  the  sample  was 
approximately  15  mW.  The  imaging  chamber  was  heated  to  maintain  the 
cells  at  37°C.  A  series  of  fluorescence  images  were  obtained  as  a  function  of 
depth  (in  0.5  pm  increments)  such  that  the  image  with  the  largest  cross- 
sectional  area  per  cell  could  be  selected  for  analysis.  At  each  depth  incre¬ 
ment,  the  cells  were  scanned  over  an  approximately  200  X  200  pm  area 
such  that  several  cells  could  be  imaged  at  once.  The  frame  rate  was  ap¬ 
proximately  1  s/slice.  The  fluorescence  images  of  the  cells  were  background 
subtracted,  Le.  by  subtracting  the  intensity  observed  in  a  region  of  the  same 
sample  where  no  cells  were  present. 

Viability  after  each  imaging  experiment  was  verified  by  checking  the 
imaged  cells  for  normal  growth  the  next  day.  The  viability  of  the  imaged 
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Table  1.  Cell  viability  data  compiled  by  cell  type,  collected  before  and 
after  each  fluorescence  experiment  for  four  of  the  five  cell  types  evaluated 
in  this  study.  The  numbers  in  parentheses  within  the  first  column 
correspond  to  the  number  of  independent  cultures  evaluated  for  each  cell 
type.  The  time  period  between  pre-  and  postexperiment  viability 
measurements  was  approximately  45  min 


Cell  line 

Viability  before 
experiments  (%) 

Viability  after 
experiments  (%) 

MCFIO  pZIP  (n  =  8) 

91.9  ±  2.5 

87.5  ±  5.3 

MCFIO  R-Ras  (n  =  7) 

93.3  ±  2.1 

84.4  ±11.3 

T47D  (n  ===  7) 

85.2  ±  1.7 

80.7  ±  4.5 

MDA231  (n  =  7) 

88.0  ±  4.8 

85.4  ±  4.7 

cells  was  compared  with  that  of  nonimaged  cells  grown  and  maintained 
under  identical  experimental  conditions.  The  imaged  cells  appeared  to 
remain  viable,  with  no  apparent  effects  caused  by  the  imaging  process.  Note 
that  the  assessment  of  viability  was  not  carried  out  quantitatively  because 
the  cells  would  have  had  to  be  detached  and  counted  in  a  hemacytometer. 
This  procedure  would  require  the  mixing  of  a  large  number  of  nonimaged 
cells  with  a  small  number  of  imaged  cells. 

RESULTS 

Table  1  shows  the  cell  viability  data  compiled  by  cell  type, 
collected  before  and  after  each  fluorescence  experiment,  for  four  of 
the  five  cell  types  evaluated  in  this  study.  The  numbers  in  paren¬ 
theses  within  the  first  column  correspond  to  the  number  of 
independent  cultures  evaluated  for  each  cell  type.  The  time  period 
between  pre-  and  postexperiment  viability  measurements  was 
approximately  45  min.  The  viability  of  all  cell  types  was  greater 
than  85%  before  the  fluorescence  experiment  and  decreased  by  less 
than  10%  during  the  course  of  the  experiment. 

Figure  1  shows  representative  fluorescence  EEM  for  the  (1) 
MCFIO  pZIP  cells  and  (2)  MCFIO  R-Ras  cells  measured  at  con¬ 
centrations  of  1.07  X  10^  cells/mL  and  1.20  X  10^  cells/mL,  re¬ 
spectively.  The  plots  are  shown  on  a  log  contour  scale,  where  each 
contour  corresponds  to  levels  of  equal  fluorescence  intensity. 
Fluorescence  peaks  “T,”  “N”  and  “F”  are  attributed  to  tryptophan, 
NAD(P)H  and  FAD,  respectively.  The  same  three  peaks  were 
observed  in  the  EEM  of  the  other  cell  types.  Thus,  all  subsequent 
analyses  focused  on  the  fluorescence  peaks  corresponding  to  these 
three  fluorophores.  However,  an  excitation-emission  wavelength 
pair  of  300-340  nm  was  used  to  quantify  tryptophan  fluorescence 
because  the  fluorescence  intensities  at  an  excitation-emission 
wavelength  pair  of  280-340  nm  (corresponding  to  the  peak  fluo¬ 
rescence  of  tryptophan)  were  saturated  at  some  of  the  higher  cell 
concentrations.  An  evaluation  of  the  fluorescence  emission  spectra 
of  tryptophan,  NAD(P)H  and  FAD  at  excitation  maxima  of  300, 
340  and  460  nm,  respectively,  indicated  that  the  spectral  line 
shapes  of  these  fluorophores  are  very  similar  for  all  cell  types. 

Figure  2  displays  the  fluorescence  intensity  vs  cell  concentration 
for  MCFIO  pZIP  and  MCFIO  R-Ras  ceUs  at  an  excitation- 
emission  maximum  of  (Fig.  2a, b)  300-340  nm  (which  corresponds 
to  the  fluorescence  of  tryptophan)  and  (Fig.  2c,d)  340-450  nm 
(which  conresponds  to  the  fluorescence  of  NAD(P)H).  The  stan¬ 
dard  deviations  of  the  fluorescence  intensity  and  the  cell  con¬ 
centration  are  shown  for  each  data  point.  A  least  squares  fit  over 
the  linear  range  is  shown  in  each  plot. 

At  MCFIO  pZIP  cell  concentrations  less  than  2.0  X  10^  cells/ 
mL,  there  is  a  linear  relationship  between  fluorescence  intensity 
and  cell  concentration  at  both  excitation-emission  wavelength 
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Figure  1.  Representative  fluorescence  EEM  for  the  MCFIO  pZIP  cells  (a) 
and  MCFIO  R-Ras  cells  (b),  measured  at  concentrations  of  1.07  X  10*^  cells/ 
mL  and  1.20  X  10^  cells/mL,  respectively.  The  plots  are  shown  on  a  log 
contour  scale,  where  each  contour  corresponds  to  levels  of  equal  fluo¬ 
rescence  intensity.  Fluorescence  peaks  “T,”  “N”  and  “F”  are  attributed  to 
tryptophan,  NAD(P)H  and  FAD,  respectively. 

pairs  given  (Fig.  2a, c).  However,  at  concentrations  that  are  greater 
than  or  equal  to  2,0  X  10^  cells/mL,  the  fluorescence  intensity 
saturates  with  increasing  cell  concentration  at  both  excitation- 
emission  wavelength  pairs  given.  This  is  presumably  because  of 
the  increased  turbidity  of  the  cell  suspension.  On  the  other  hand, 
the  fluorescence  intensity  of  MCFIO  R-Ras  cells  increases  linearly 
with  concentration  over  the  entire  range  up  to  3  X  10^  cells/mL  at 
both  excitation-emission  wavelength  pairs  (Fig.  2b,d). 

Table  2  shows  the  cell  concentration  at  which  the  linear 
relationship  between  the  fluorescence  intensity  and  ceU  concentra¬ 
tion  saturates  for  four  of  the  five  cell  types  evaluated  in  this  study. 
These  maximum  cell  concentration  values  were  found  to  be  the 
same  for  fluorescence  measurements  at  the  three  excitation-emis¬ 
sion  maxima  corresponding  to  tryptophan,  NAD(P)H  and  FAD. 
The  normal  MCFIO  pZIP  cell  line  displays  saturation  in  fluo¬ 
rescence  intensity  at  a  lower  concentration  than  that  of  the  three 
malignant  cell  types. 

Figure  3  shows  the  slope  and  the  standard  error  for  four  of  the 
five  cell  types  evaluated  in  this  study  at  the  three  excitation- 
emission  maxima  corresponding  to  the  fluorescence  of  tryptophan. 
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Figure  2.  Fluorescence  intensity  vj  cell 
concentration  for  MCFIO  pZIP  and 
MCFIO  R-Ras  cells  at  an  excitation- 
emission  maximum  of  300-340  nm 
(which  corresponds  to  the  fluorescence 
of  tryptophan)  (a  and  b,  respectively)  and 
34(M50  nm  (which  corresponds  to  the 
fluorescence  of  NAD(P)H)  (c  and  d, 
respectively).  The  standard  deviations  of 
the  fluorescence  intensity  and  cell  con¬ 
centration  are  shown  for  each  data 
point.  A  least  squares  fit  over  the  linear 
range  is  shown  in  each  plot. 
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NAD(P)H  and  FAD.  The  slope  corresponds  to  the  fluorescence 
intensity  per  cell.  The  slopes  are  multiplied  by  a  scaling  factor  to 
display  the  fluorescence  per  cell  of  the  three  fluorophores  on  the 
same  graph  (shown  in  legend). 

General  linear  hypothesis  testing  (28)  was  used  to  determine  the 
level  of  significance  at  which  differences  in  tryptophan,  NAD(P)H 
and  FAD  fluorescence  are  observed  between  the  normal  and  malig¬ 
nant  cells.  All  the  malignant  cells  showed  a  statistically  significant 
decrease  in  the  tryptophan  fluorescence  per  cell  relative  to  that  of 
the  normal  MCFIO  pZIP  cells,  at  a  level  of  significance  of  F  < 
0.002.  No  statistically  significant  differences  were  observed  in  the 
NAD(P)H  fluorescence  per  cell  between  the  normal  MCFIO  pZIP 
cells  and  any  of  the  malignant  cell  types.  The  FAD  fluorescence 
per  cell  of  the  MCFIO  R-Ras  cells  was  increased  relative  to  that  of 
the  MCFIO  pZDP  cells  at  a  level  of  significance  of  P  <  0.05. 
However,  no  statistically  significant  differences  were  observed  in 
the  FAD  fluorescence  per  cell  between  the  normal  MCFIO  pZIP 
cells  and  any  of  the  genetically  dissimilar  malignant  cell  types. 

As  mentioned  earlier  in  Materials  and  Methods,  the  fluorescence 
intensity  vs  cell  concentration  was  also  measured  from  the 
malignant  T47D  cells  transfected  with  the  R-Ras  oncogene.  In 
this  case,  tryptophan  and  NAD(P)H  fluorescence  per  cell  was  not 
significantly  altered,  but  a  statistically  significant  increase  in  FAD 
fluorescence  was  observed  in  the  T47D  R-Ras  cells  relative  to  that 
in  the  T47D  cells.  These  findings  indicate  that  the  decrease  in 
tryptophan  fluorescence  observed  with  transformation  is  not 
attributed  to  the  expression  of  the  R-Ras  oncogene  alone  but 
rather  to  the  transformation  of  the  cell  from  a  normal  to  a  malignant 
phenotype.  However,  the  increase  in  FAD  fluorescence  in  the 
MCFIO  R-Ras  cells,  relative  to  that  in  the  MCFIO  pZIP  cells,  may 


be  because  of  the  transfection  of  the  R-Ras  oncogene  rather  than 
due  to  transformation  from  a  normal  to  a  malignant  phenotype. 

Figure  4  shows  the  average  cell  cross-sectional  areas  (de¬ 
termined  from  approximately  20  cells  for  each  cell  type)  for  four  of 
the  five  cell  types  evaluated  in  this  study.  It  can  be  seen  that  all  of 
the  cell  types  are  similar  in  size,  with  cell  sizes  ranging  from  150  to 
200  pm^.  The  variability  in  size  within  a  cell  type  appears  to  be 
typical  of  these  cells.  Figure  4  indicates  that  the  differences  in 
fluorescence  intensity  per  cell  between  the  normal  and  malignant 
cells  are  not  attributed  to  differences  in  cell  size. 

Finally,  NAD(P)H  fluorescence  obtained  from  cells  in  suspen¬ 
sion  were  compared  with  that  obtained  from  cell  monolayers 
imaged  using  two-photon  fluorescence  microscopy.  Figure  5  shows 
representative  fluorescence  images  of  (1)  MCFIO  pZIP  and  (2) 
MCFIO  R-Ras  cell  monolayers  at  775  nm  excitation  (corresponds 
to  single-photon  excitation  of  387.5  nm),  imaged  under  identical 


Table  2.  The  cell  concentration  at  which  the  linear  relationship  between 
the  fluorescence  intensity  and  the  cell  concentration  saturates  for  four  of  the 
five  cell  types  evaluated  in  this  study.  These  cell  concentration  values  were 
found  to  be  the  same  for  fluorescence  measurements  at  the  three  excitation- 
emission  maxima  corresponding  to  tryptophan,  NAD(P)H  and  FAD 


Cell  line 

Cell  concentration  at  which 
fluorescence  saturates  (cells/mL) 

MCFIO  pZIP 

2.0  X  10® 

MCFIO  R-Ras 

>3  X  10® 

T47D 

>3  X  10® 

MDA231 

>3  X  10® 
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Figure  3.  Slope  and  standard  error  for  four  of  the  five  cell  types  evaluated 
in  this  study  at  the  three  excitation-emission  maxima  corresponding  to  the 
fluorescence  of  tryptophan,  NAD(P)H  and  FAD.  The  slope  corresponds  to 
the  fluorescence  intensity  per  cell.  The  slopes  are  multiplied  by  a  scaling 
factor  to  display  the  fluorescence  per  cell  of  the  three  fluorophores  on  the 
same  graph  (shown  in  legend). 


experimental  conditions.  It  can  be  seen  that  the  fluorescence  of  the 
MCFIO  R-Ras  cells  is  approximately  a  factor  of  2  lower  than  that 
of  the  MCFIO  pZIP  cells.  To  evaluate  these  results  quantitatively, 
Image!  (NIH,  http://rsb.info.nih.gov/ij)  was  used  to  (1)  manually 
identify  the  fluorescence  image  that  displayed  the  largest  cross- 
sectional  area  per  cell;  (2)  integrate  the  fluorescence  intensity  from 
each  cell;  and  (3)  average  the  fluorescence  intensity  from  a  total 
of  five  cells  within  each  monolayer.  The  ratio  of  the  average 
fluorescence  intensity  of  the  MCFIO  pZIP  and  MCFIO  R-Ras  cell 
monolayers  (imaged  under  identical  experimental  conditions)  was 
then  calculated.  This  ratio  was  averaged  for  three  independent  pairs 
of  MCFIO  pZIP  and  MCFIO  R-Ras  monolayer  images  and 
calculated  to  be  2.3  ±  1.0.  A  paired  Mest  (27)  was  used  to  deter¬ 
mine  the  level  of  significance  at  which  differences  in  NAD(P)H 
fluorescence  are  observed  between  the  MCFIO  pZIP  and  MCFIO 
R-Ras  monolayers.  The  NAD(P)H  fluorescence  of  the  MCFIO  R- 
Ras  cells  was  decreased  relative  to  that  of  the  MCFIO  pZIP  cells  at 
a  level  of  significance  of  P  <  0.05.  In  addition,  a  qualitative  ex¬ 
amination  of  the  images  showed  that  there  are  differences  in  the 
spatial  distribution  of  NAD(P)H  fluorescence  between  the  MCFIO 
pZIP  and  MCFIO  R-Ras  ceUs.  Specifically,  the  NAD(P)H  fluo¬ 
rescence  of  the  normal  MCFIO  pZIP  cells  appeared  to  be  more 
granular  than  that  of  the  transformed  MCFIO  R-Ras  cells. 

DISCUSSION 

The  ceU  suspension  studies  show  that  there  is  a  statistically 
significant  decrease  in  the  tryptophan  fluorescence  per  ceU  in  the 
collective  set  of  malignant  cell  types  compared  with  that  of  the 
normal  cell  line.  However,  no  statistically  significant  differences 
were  observed  in  the  NAD(P)H  or  FAD  fluorescence  per  cell 
between  the  normal  cells  and  any  of  the  malignant  cell  types,  based 
on  the  cell  suspension  studies.  Two-photon  fluorescence  micros¬ 
copy  of  cell  monolayers  indicated  a  statistically  significant  de¬ 
crease  in  the  NAD(P)H  fluorescence  of  the  MCFIO  R-Ras  cells 
relative  to  that  of  the  MCFIO  pZIP  cells  (normal  and  malignant 
cells  of  similar  genetic  origins).  These  results  suggest  that  fluo¬ 
rescence  imaging  of  single  cells  or  the  cell  monolayer  preparation 
may  provide  more  contrast  than  volume-averaged  fluorescence 
measurements  of  cells  in  suspension.  The  ultimate  goal  of  this 
work  is  to  apply  these  findings  to  observations  in  tissue.  However, 
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Figure  4.  Average  cell  cross-sectional  areas  (determined  from  approxi¬ 
mately  20  cells  for  each  cell  type)  for  four  of  the  five  cell  types  evaluated  in 
this  study. 


simple  cell  studies  such  as  the  one  described  here  are  not  directly 
comparable  with  tissue  studies  because  a  number  of  other  factors 
are  present  in  tissue,  such  as  the  interaction  with  the  extracellular 
matrix  and  vasculature.  Therefore,  it  is  premature  to  make  any 
links  between  the  results  of  our  ceU  studies  and  previous  ex  vivo 
studies  of  the  breast. 

The  differences  in  the  intrinsic  fluorescence  properties  of  the 
normal  and  malignant  cells  reflect  the  amount  of  fluorophore 
present  or  the  properties  of  the  surrounding  microenvironment  (or 
both)  (7).  The  decreased  tryptophan  fluorescence  observed  in  the 
malignant  cell  lines  may  indicate  that  there  is  a  lower  protein 
content  in  the  malignant  cells  relative  to  that  of  the  normal  cells.  In 
addition,  a  decreased  quantum  yield  could  result  from  changes  in 
the  microenvironment,  such  as  a  change  in  the  local  prevalence  of 
quenchers  such  as  oxygen. 

The  results  reported  for  NAD(P)H  in  our  study  are  inconsistent 
with  that  reported  by  Glassman  et  al.  (22),  who  found  that  there  is 
an  increase  in  NAD(P)H  fluorescence  with  malignancy.  In  our  cell 
suspension  measurements,  no  statistically  significant  differences  in 
the  NAD(P)H  fluorescence  was  observed  between  the  normal  and 
any  of  the  malignant  cells,  whereas  in  our  two-photon  fluorescence 
microscopy  study,  there  was  a  statistically  significant  decrease  in 
the  NAD(P)H  fluorescence  of  the  malignant  cells,  compared  with 
their  genetically  similar  normal  counterpart.  The  inconsistency  be¬ 
tween  our  study  and  that  reported  by  Glassman  et  al,  (22)  could  be 
attributed  in  part  to  differences  in  the  study  design  and  the  sample 
preparation  protocol. 

It  is  instructive  to  compare  the  findings  from  our  cell  suspension 
and  cell  monolayer  studies  with  those  that  used  similar  protocols 
(evaluation  of  normal  and  malignant  cell  lines  of  genetically 
similar  origins,  viability  testing  and  fluorescence  measurements  at 
nonturbid  cell  concentrations).  Grossman  et  al.  (18)  characterized 
the  fluorescence  properties  of  normal  and  malignant  murine 
fibroblast  cell  suspensions;  the  malignant  cell  line  was  obtained 
by  transforming  the  normal  cell  line  with  the  H-Ras  oncogene. 
They  observed  a  decrease  in  tryptophan  and  NAD(P)H  fluores¬ 
cence  with  malignancy.  Their  observations  relating  to  tryptophan 
fluorescence  are  consistent  with  that  observed  in  our  cell  sus¬ 
pension  studies,  and  their  findings  related  to  NAD(P)H  fluores¬ 
cence  are  consistent  with  that  observed  in  our  cell  monolayer 
studies.  Pitts  et  al.  (5)  examined  the  fluorescence  properties  of 
human  bronchial  epithelial  cells,  including  a  normal  cell  line  and 
a  tobacco  carcinogen-transformed  cell  line  of  the  same  genetic 
origin.  They  found  that  the  NAD(P)H  and  FAD  fluorescence 
decreased  with  malignancy  but  that  the  tryptophan  fluorescence 
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remained  unchanged.  Their  results  for  NAD(P)H,  which  they 
obtained  from  cell  monolayers,  are  consistent  with  those  observed 
in  our  cell  monolayer  studies.  Their  results  for  FAD  and  trypto¬ 
phan  are  inconsistent  with  our  findings.  These  inconsistencies  may 
be  attributed  to  the  method  of  transformation  of  the  cells,  the 
cell  lines  used,  the  cell  preparation  method  or  the  measurement 
technique. 

Our  results  are  inconsistent  with  those  reported  from  a  recent 
study  on  cervical  tissue  slices  (29),  which  showed  that  there  is  an 
increase  in  NAD(P)H  fluorescence  with  neoplasia.  The  inconsis¬ 
tency  between  the  tissue  slice  study  and  our  cell  study  may  be 
explained  by  the  fact  that  the  NAD(P)H  fluorescence  will  be 
influenced  not  only  by  the  metabolic  demand  but  also  by  the 
oxygen  supply  available  to  the  cells.  The  metabolic  demand  is 
higher  for  the  proliferative  malignant  cells  compared  with  their 
normal  counterpart  in  both  cells  (30)  and  tissues  (31,32).  However, 
the  oxygen  supply  may  be  quite  different  in  the  two  systems.  The 
cellular  microenvironment  is  hyperoxic  compared  with  the  tissue 
microenvironment.  Thus,  whereas  an  increase  in  the  metabolic 
demand  may  increase  NAD(P)H  concentration  (and  fluorescence) 
in  tissues,  the  reverse  may  be  true  in  cells,  where  the  oxygen  supply 
is  not  constrained.  As  will  be  discussed  below,  more  complex  bio¬ 
logical  models  can  be  used  in  the  future  to  evaluate  this  effect. 

The  normal  MCFIO  pZIP  cell  line  showed  saturation  in 
fluorescence  intensity  at  concentrations  greater  than  or  equal  to 
2.0  X  10^  cells/mL,  whereas  each  of  the  malignant  cell  lines  showed 
a  linear  relationship  between  fluorescence  intensity  and  cell 
concentrations  up  to  3  X  10^  cells/mL.  The  saturation  in  fluo¬ 
rescence  intensity  is  likely  caused  by  the  increased  turbidity  of  the 
cell  suspension.  Therefore,  these  results  suggest  that  the  MCFIO 
pZIP  cell  suspension  is  more  turbid  than  the  malignant  cell 
suspensions  at  similar  cell  concentrations.  This  may  indicate  that 
either  the  absoiption  or  the  scattering  properties  of  normal  cells  are 
increased  relative  to  those  of  the  malignant  cells.  Because  saturation 
of  the  fluorescence  intensity  occurs  at  the  same  MCFIO  pZIP  cell 
concentration  at  all  three  excitation-emission  maxima,  the  non¬ 
linearity  in  fluorescence  intensity  concentration  is  more  likely 
because  of  scattering  (rather  than  due  to  absorption).  Therefore,  it  is 
hypothesized  that  the  MCFIO  pZIP  cells  scatter  more  than  their 
malignant  cell  counterpart.  There  are  several  studies  that  support 
this  hypothesis.  Croce  et  al.  (33)  have  previously  reported  that 
fibroblasts  with  a  malignant  phenotype  exhibit  a  decrease  in  the 
number  of  mitochondria  and  a  loss  of  organization  relative  to 
a  putatively  normal  ceU  line  (this  may  also  explain  the  decreased 
granularity  seen  in  the  mitochondrial  NAD(P)H  fluorescence  in  the 
transformed  cells  in  Fig.  5).  Because  the  mitochondria  contribute  to 
a  substantial  amount  of  cellular  scattering  (34),  the  malignant  cells 
will  have  decreased  scattering  relative  to  their  normal  counterpart. 
Zonios  et  al.  (35)  observed  a  decrease  in  the  scattering  coefficient  in 
adenomatous  colon  polyps  in  vivo,  relative  to  normal  tissues,  at 
wavelengths  less  than  550  nm.  Georgakoudi  et  al.  (36)  also  found 
that  scattering  decreased  with  disease  progression  in  patients  with 
Barrett  esophagus. 

Backman  et  al.  (37)  report  that  there  is  increased  scattering  from 
neoplastic  cells  relative  to  nonneoplastic  cells  and  show  that  this  is 
due  to  increased  scattering  from  the  enlarged  nucleus  of  the  neo¬ 
plastic  cells.  Their  method  examined  polarized  backscattered  light, 
which  they  report  is  primarily  due  to  nuclear  size  (they  report  that 
nuclear  scattering  accounts  for  the  majority  of  small-angle  forward 
scattering  and  backscattering,  whereas  high-angle  scattering  is  due 
to  smaller  organelles  such  as  the  mitochondria)  (37).  In  our  study, 
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Figure  5.  Representative  fluorescence  image  of  MCFIO  pZIP  (a)  and 
MCFIO  R-Ras  (b)  cell  monolayers  at  775  nm  excitation  (corresponds  to 
single-photon  excitation  of  387.5  nm),  imaged  under  identical  experimental 
conditions.  Please  note  the  different  scales  in  (a)  and  (b). 

the  transformed  and  nontransformed  cells  had  similar  nuclear  sizes; 
thus,  differences  in  scattering  may  be  primarily  attributed  to  the 
differences  in  the  mitochondrial  content  of  malignant  and  non- 
malignant  breast  cells. 

In  summary,  significant  differences  have  been  noted  in  the  fluo¬ 
rescence  properties  of  normal  and  malignant  human  breast  cells.  In 
future  studies  it  may  be  possible  to  build  upon  these  findings  in 
order  to  understand  how  the  contribution  of  cellular  fluorescence 
affects  fluorescence  spectra  measured  in  intact  tissue.  For  example, 
it  is  possible  to  grow  epithelial  cells  in  a  three-dimensional  colla¬ 
gen  matrix,  thus  producing  a  system  in  which  cell  growth  and 
differentiation  are  modulated  by  the  effects  of  the  extracellular 
matrix  (38).  Using  such  a  system,  one  could  potentially  examine 
cellular  fluorescence  with  the  added  influence  of  the  extracellular 
matrix  on  cell  growth  and  differentiation.  To  more  accurately 
simulate  in  situ  conditions,  normal  and  malignant  human  breast 
cell  lines  can  be  grown  subcutaneously  in  a  nude  mouse  model 
(39).  In  such  a  system  the  added  effects  of  the  vasculature,  in 
addition  to  the  extracellular  matrix,  on  the  tissue  fluorescence  could 
be  studied.  A  stepwise  approach  such  as  this  presents  the  possi¬ 
bility  to  bridge  these  studies  to  fluorescence  measurements  made 
on  intact  human  breast  tissue. 
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Abstract — Nonmalignant  (n  =  36)  and  malignant  (n  =  20) 
tissue  samples  were  obtained  from  breast  cancer  and  breast 
reduction  surgeries.  These  tissues  were  characterized  using  mul¬ 
tiple  excitation  wavelength  fluorescence  spectroscopy  and  diffuse 
reflectance  spectroscopy  in  the  ultraviolet-visible  wavelength 
range,  immediately  after  excision.  Spectra  were  then  analyzed 
using  principal  component  analysis  (PCA)  as  a  data  reduction 
technique.  PCA  was  performed  on  each  fluorescence  spectrum,  as 
well  as  on  the  diffuse  reflectance  spectrum  individually,  to  estab¬ 
lish  a  set  of  principal  components  for  each  spectrum.  A  Wilcoxon 
rank-sum  test  was  used  to  determine  which  principal  components 
show  statistically  significant  differences  between  malignant  and 
nonmalignant  tissues.  Finally,  a  support  vector  machine  (SVM) 
algorithm  was  utilized  to  classify  the  samples  based  on  the 
diagnostically  useful  principal  components.  Cross-validation  of 
this  nonparametric  algorithm  was  carried  out  to  determine  its 
classification  accuracy  in  an  unbiased  manner.  Multiexcitatiou  f 
fluorescence  spectroscopy  was  successful  in  discriminating  m^g* 
nant  and  nonmalignant  tissues,  with  a  sensitivity  and  specifidt^  of 
70%  and  92%,  respectively.  The  sensitivity  (30%)  and  speclfldty 
(78%)  of  diffuse  reflectance  spectroscopy  alone  was  significantly 
lower.  Combining  fluorescence  and  diffuse  reflectance^pectra  did 
not  improve  the  classification  accuracy  of  an  algoritto  based  on 
fluorescence  spectra  alone.  The  fluorescence  excit^flon-emission 
wavelengths  identified  as  being  diagnostic  from 4he  PCA-SVM  al¬ 
gorithm  suggest  that  the  important  fluoropho]|||for  breast  <^cer 
diagnosis  are  most  likely  tryptophan,  NAD(P)ll^id  ^^^Ot^^bteins. 

Index  Terms — Breast  cancer,  fluorescence,  reflecMce,  spec¬ 
troscopy. 

L  Introduction 

Breast  cancer  screening  relies  on  mamrdbgraphy  and 
clinical  breast  exam  for  identifying  suspicious  lesions. 
Manunograms  are  useful  in  detecting  small  tumors  that  may  4 
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be  unnoticeable  by  physical  exam;  however,  they  produce 
a  significant  percentage  of  false  positives  [1].  As  a  result, 
follow  up  diagnostic  procedures  such  as  percutaneous  needle 
biopsy  are  performed  to  further  evaluate  breast  abnormalities 
[1],  Needle  biopsy  has  a  limited  sampling  accuracy  because 
only  a  few  small  pieces  of  tissue  are  extracted  from  random 
locations  in||ie  suspicious  mass.  In  some  cases,  sampling  of 
the  suspicious  mass  may  be  missed  altogether.  Consequences 
include  ^false-negative  rate  of  \%-l%  [2]  when  verified  with 
followup  mammography,  as  well  as  the  requirement  of  repeat 
biopsies  (percut&eous  or  surgical)  in  9%-18%  of  patients 
[3],  [4]  (due  to; discordance  between  histplogical  findings  and 
fmammdpa^y).  jf /  ' 


core^jpedle  enabling  sampling  O^ift^tiple  tissue  sites 

witnout  the  need  for  biopsy,  and  by  providinS|f  al-time  diag¬ 
nosis.  Fluoresbience  and  diffuse  reflectsmce  specWscopic  mea¬ 
surements  Mthe  UV-VIS  have  beeh  showhtb^have  consistently 
high  sensiflvities  and  specificitieg  for  detectMf  buman  pre-can¬ 
cers  anil  cancers,  in  vivo  in  a  ni|rter  of  orgaj||sites  including 
thegisffointestinal  tract,  bronchulpind  cervix^.  Additionally, 
several  groups  have  demofisK^d  that tber^ire  significant  dif¬ 
ferences  in  the  fluores&nce  anidif^se  reflectance  spectra  of 
normal,  benign  and  nfSignant  br^SJt  tissues  [6]-[16].  Incorpo¬ 
rating  UV-VIS  specti^copy  as  an|^unct  diagnostic  modality 
to  core  needle  biopsy  %puld  potentially  reduce  the  false  nega¬ 
tive  rate  of  cuTOht  needle  tiiopsyjprocedures.  This  in  turn  could 


potefitMly  lead-tQ  fewer  repeat  biopsies  in  patients  with  suspi- 
ious  breast  lesidns. 

The  end^|enous  fluorophores  present  in  breast  tissue  include 
tryptophai^leduced  nicotinamide  adenine  dinucleotide  (phos¬ 
phate)  p|Ab(P)H),  flavin  adenine  dinucleotide  (FAD),  and  col- 
Tryptophan  fluorescence,  which  has  an  excitation, 
enussibn  maximum  of  300, 340  nm,  is  an  indicator  of  protein  or 
free  t^tophan  content.  NAD(P)H  and  flavoproteins  are  indi¬ 
cators  of  metabolic  activity  [17]  and  have  excitation,  emission 
maxima  at  35 1 , 460  nm  and  450-520  nm,  respectively  [17].  Col¬ 
lagen  is  the  primary  structural  protein  in  the  extracellular  matrix. 
It  has  several  excitation,  emission  maxima,  one  of  which  occurs 
at  325,  390  nm  (due  to  cross-link  fluorescence)  [17]. 

Alfano  et  al  [12]  were  the  first  to  measure  fluorescence 
spectra  of  malignant  and  nonmalignant  breast  tissues.  Since 
then,  several  groups  have  investigated  the  use  of  fluorescence 
spectroscopy  for  breast  cancer  detection  in  ex  vivo  studies 
[6]-[12].  In  the  collective  set  of  studies  to  date,  emission 
spectra  were  acquired  at  excitation  wavelengths  of  300,  340, 
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458, 488,  and  514  nm  and  excitation  spectra  were  measured  at 
emission  wavelengths  of  340,  390,  and  460  nm  from,  normal, 
benign,  and  malignant  breast  tissues.  Emission  spectroscopy 
measures  the  fluorescence  spectrum,  while  excitation  spec¬ 
troscopy  measures  the  absorption  spectrum  of  a  fluorophore. 
The  wavelength-dependent  effects  of  nonfluorescent  absorbers 
and  scatterers  in  tissue  also  affect  the  tissue  excitation  and 
emission  spectra, 

Gupta  [7]  and  Majumder  et  al  [6]  analyzed  different  data  sets 
collected  from  the  same  set  of  breast  tissues  ex  vivo  and  showed 


cence  and  diffuse  reflectance  spectra  in  order  to  classify  a  partic¬ 
ular  sample  as  malignant  or  nonmalignant  (benign,  normal)  and 
a  cross-validation  scheme  was  utilized  to  test  the  algorithm’s 
classification  accuracy  in  an  unbiased  manner. 

n.  Methods 

A.  Tissue  Handling  Protocol  for  Ex  Vivo  Studies 

All  of  the  fluorescence  and  most  of  the  diffuse  reflectance 
measurements  to  date  have  been  made  on  excised  breast  tissues. 


that  the  emission  spectra  at  excitation  wavelengths  of  340  and 
488  nm  and  excitation  spectra  at  emission  wavelengths  390  and 
460  nm  exhibit  significant  differences  between  normal,  benign 
and  malignant  tissues.  Using  the  integrated  emission  intensities 
at  an  excitation  wavelength  of  340  nm,  in  a  binary  classification 
scheme,  they  were  able  to  differentiate  malignant  from  normal 
and  benign  tissues  with  a  sensitivity  and  specificity  of  98%. 
The  fluorescence  was  attributed  to  reduced  nicotinamide  ade¬ 
nine  dinucleotide  (NADH)  and  collagen.  Yang  et  al  [8]-[10] 
reported  that  emission  spectra  at  300-nm  excitation  and  exci¬ 
tation  spectra  at  340-nm  emission  (tryptophan  is  the  primary 
fluorophore)  can  separate  malignant  and  nonmalignant  tissues. 
They  found,  for  example,  that  the  ratio  of  normalized  intensi¬ 
ties  at  289  and  268  nm  of  the  excitation  spectra  discriminate  i 
between  malignant  and  fibrous  tissues  with  93%  sensitivity  and® 
95%  specificity.  Although  previous  work  shows  promising;  fe- 
suits,  there  are  important  gaps  in  these  studies  that  still  r^d  to 
be  addressed.  The  main  limitations  of  previous  studies  ^®IT) 
fluorescence  spectra  were  obtained  only  at  one  or  several  exMa® . 
tion  wavelengths  from  breast  tissues;  2)  the  utility  of?e6mbinidg^|| 
fluorescence  and  diffuse  reflectance  spectroscop^lwas  not  eval¬ 
uated;  and  3)  the  classification  accuracy  of  the  Ifiagnostic  „^go- 
rithms  were  not  tested  in  an  unbiased  manner. 

Spectral  differences  observed  in  the  fluorescencd^pectra 
of  normal,  benign  and  malignant  br^st  tissuefe^^^  also  be  j 
attributed  in  part  to  nonfluorescent  abs^^bers  a^iscatterersA!' 
Diffuse  reflectance  spectroscopy  providH^direct  ni'eayemefetr 
of  tissue  absorption  as  well  as  scattering.  Tfe  prima^^onllu- 
orescent  absorbers  in  the  UV-VIS  in  breaSi^ssue  are^DNA, 
oxygenated  and  deoxygenated  hemoglobin  and  /3-carotene 

Several  groups  have  explored  the  utility  of  diffuse  reflectance 
spectroscopy  between  250-800  nm,  for  breast  cancer  detection 
ex  vivo  [10],  [13H15]  and  in  v/vo[i5],  [16].  Bigio  et  al  [I6f  l 
measured  diffuse  reflectance  spectra  through  a  core  biopsy 
needle  and  during  breast  cancer  surgery  and  showed  that  this 
technique  can  differentiate  malignant  from  normal  tissues  with 
a  sensitivity  of  60%-70%  and  a  specificity  of  85%-95%.  The 
collective  set  of  ex  vivo  and  in  vivo  studies  [10],  [13]-[16]  show 
that  diffuse  reflectance  spectra  can  differentiate  malignant 
from  normal  tissues.  The  changes  in  the  diffuse  reflectance  of 
malignant  tissues,  was  attributed  to  increased  DNA,  protein 
and  hemoglobin  absorption  and  increased  scattering. 


However,  no  systematic  investigations  have  been  carried  out  to 
assess  the  effect  of  excision  and  time  after  excision  on  the  tissue 
fluorescence  and  diffuse  reflectance  spectra.  Therefore,  the  ef¬ 
fects  of  tissue  handling  and  storage  on  optical  spectroscopy  of 
excised  tissue  specimens  were  systematically  evaluated  in  a  pre¬ 
vious  study  [19].  In  that  study,  multiexcitation  fluorescence  and 
diffuse  reflectance  spectra  of  the  epithelial  tissue  of  the  hamster 
cheek  poucBIwas  characterized  in  vivo,  ex  vivo  (immediately 
after  excisfSn)  and  after  being  frozen  and  thawed.  The  results 
indicatefThat  the  freezing  and  thawing  process  of  tissue  pro¬ 
duces  usignificantp^viation  in  intensity  and  line  shape  relative 
to  thfe  in  vivo  spectra  over  the  entire  UV-VIS  range.  On  the  other 
hand^i^  v/va^ctra  measured  within  a  short  time  after  excision 
'  were  found  fSprovide  a  rlfatively  closejij^roximation  of  in  vivo 
"spectra.  /■:, 


B,  Optical  SpietiS^kdpy  of  Breast  Tissudsf 

4;^Pproval  fgr  the  optical  spectroscopy  study  was  obtained 
Aom  the  Insrffiitional  Review  BoardI|4  the  Univ^sity  of  Wis- 
gconsin,  Madison.  Tissue  was  obtained  frorh  jatientk^dergoing 
l^^er  a  lUlipectomy,  mastectomy  or  breast  r&gction  surgery.  In 
t%i§^||^rie  clinical  protocol,  tis|ue  specimei^^iWere  harvested 
duSnt  surgery  and  normal jmd  tumor  tissue^mples  were  cut 
from  the  larger  specimenfflie  tissuS^mpIes  that  were  cut  from 
gthe  larger  specimen  \y^e  generaUy jt  least  1  cm  in  length/width 
kmd  3  mm  thick,  thu$ ;  providing^  ;Semi-infinite  geometry  for 
optical  spectros^Jy.  C^tical  spectiroscopic  measurements  were 
carried  out  on  tne  fresh  tissue  samples,  within  two  hours  after 
surgical  excisioigfMeasufemeiits  were  made  on  a  total  of  56  dif- 
fejrAti&ue  sah^les,  from  a  total  of  32  different  patients.  In  the 
majority  of  patieifts  who  participated  in  this  study,  a  normal  and 
^mor  pair|vas  obtained  for  optical  spectroscopy.  When  normal 
f  and  tum^issue  sites  were  available  from  a  patient,  the  normal 
'vtissue  ilicimen  was  obtained  from  the  same  breast  that  con- 
^flmed  toe  tumor  (bilateral  cancer  surgeries  are  relatively  rare). 
Inl^pSt  cases,  the  normal  tissue  was  obtained  as  far  as  possible 
frorri'^^e  tumor  site  (generally  1  cm  or  more  from  the  tumor 
boundal^)  m  order  to  minimize  the  potential  effects  of  malig¬ 
nancy.  In'%  minority  of  the  patients,  only  either  a  normal  or 
tumor  sample  was  available  for  the  measurement,  usually  due 
to  limitations  in  the  size  of  the  surgical  specimen. 

All  measurements  were  made  using  a  Skinscan  spectrofluo- 
rometer  (JY  Horiba,  Edison,  NJ).  This  instrument  consists  of 


The  primary  goal  of  the  study  described  here  was  to  charac¬ 
terize  the  muitiexcitation  fluorescence  spectra  (at  a  total  of  nine 
excitation  wavelengths  in  the  UV-VIS)  and  UV-VIS  diffuse  re¬ 
flectance  spectra  of  malignant  and  nonmalignant  breast  tissues 
and  to  identify  the  optimal  spectral  features  for  breast  cancer  di¬ 
agnosis.  A  novel  nonparametric  algorithm  (applicable  for  small 
sample  sizes)  was  developed  to  analyze  the  breast  tissue  fluores¬ 


a  150  W  xenon  lamp,  double  grating  excitation  and  emission 
monochromators  having  fixed  bandpasses  of  5  nm  for  both  ex¬ 
citation  and  emission,  and  a  photomultiplier  tube  (PMT).  The 
adjustable  parameters  of  the  system  are  the  wavelength  range 
and  increment,  the  signal  acquisition  time  and  the  PMT  high 
voltage.  The  illumination  and  collection  of  light  was  coupled 
through  a  fiber  optic  probe,  consisting  of  a  central  collection 
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core  with  a  diameter  of  1.52  mm  surrounded  by  an  illumination 
ring,  with  an  outer  diameter  of  2.18  mm.  Both  the  illumination 
ring  and  collection  core  are  made  up  of  31  individual  fibers, 
each  with  a  core/cladding  diameter  of  200/245  ^m.  The  numer¬ 
ical  aperture  of  the  excitation  and  emission  fibers  is  0.22  and 
0.12,  respectively.  The  probe  geometry  was  designed  to  max¬ 
imize  the  probing  depth  in  the  tissue,  while  providing  good 
signal-to-noise  ratio.  Simulations  using  a  modified,  three-di¬ 
mensional,  weighted-photon  Monte  Carlo  code  have  been  car- 


TABLE  I 

Histological  Breakdown  of  the  56 
Breast  Tissue  Samples  Examined  in  the  Optical  Spectroscopy  Study 


Totnl  Tissues 

20 

Invasive  ductal  cm'cinoina  (IDC) 

16 

Invasive  lobuiai'  cai'cinoina 

Carcinoma  in  situ  (CIS) 

1 

CIS  +  IDC  combined 

1 

Norma FBeiiigtt  Fibrous  Tissues  1 5 


ried  out  to  evaluate  the  probing  depth  achieved  with  this  illu¬ 
mination  and  collection  geometry  in  a  turbid  medium  [20].  For 
a  homogeneous  medium  with  a  fixed  scattering  coefficient  of 
1 10.4  cm“^,  and  absorption  coefficient  varying  from  31.8  cm“^ 
to  1.3  cm“^  (i.e.,  absorption  coefficients  spanning  the  UV-VIS 
range),  the  probing  depth  varied  from  450  to  1350  pun.  In  a 
turbid  medium  with  a  fixed  absorption  coefficient  of  10.8  cm“^, 
and  scattering  coefficient  varying  from  225  cm“^  to  50  cm'^^ 
(i.e.,  scattering  coefficients  spanning  the  UV-VIS  range),  the 
probing  depth  varied  from  550  to  1050  pbm.  The  probing  depth 
in  a  medium  with  absorption  (2.3  cm“^)  and  scattering  coeffi¬ 
cients  (167  cm“^)  representative  of  breast  adipose  tissue  at  540 
nm  [21]  was  1050  p,m. 

Fluorescence  emission  spectra  were  recorded  at  excitation 
wavelengths  of  300  to  460  nm,  in  20-nm  increments.  At  each  ex-  , 
citation  wavelength,  fluorescence  emission  was  recorded  from; 
10  nm  greater  than  the  excitation  wavelength,  up  to  600  nn^^ 
in  5-nm  increments  (e.g.,  310  to  600  nm  at  an  excitation  Waye- 
length  of  300  nm).  Fluorescence  emission  spectra  were  obra&ed 
at  a  total  of  nine  excitation  wavelengths.  A  signal  aqquisitio^^ 
time  of  0.5  s  was  used  per  excitation-emission  wavelength  pair 
and  the  PMT  high  voltage  was  set  at  950  V  for  all  measurements. 
The  total  measurement  time  was  approximately;^  min.  The  flu¬ 
orescence  emission  spectral  intensities  at  eaph  excitation;  Wave¬ 
length  were  normalized  to  that  of  a  reference^jpftptodiod^  which 
accounts  for  wavelength-dependent  variations  the  Excitation 

light  intensity.  The  fluorescence  emissmli  spectniimat  each  ex- i 
citation  wavelength  was  corrected  for  wavelength  dependencig^ 
in  the  emission  monochromator  and  PMT  using  instrument  SfE- 
cific  correction  factors.  Finally,  all  fluorescence  emission  spec¬ 
tral  intensities  were  divided  by  the  peak  fluorescence  inten¬ 
sity  (excitation-emission:  460-575  nm)  of  a  Rhoda^pe  B  dye 
(1 15H3423,  Sigma  Chemical  Co.,  MO)  dissolved  in  Hater  (2. 14 
/iM),  measured  after  each  patient  study  to  correct  for  time-de/ 
pendent  changes  in  the  throughput  of  the  instrument.  " 

Diffuse  reflectance  was  measured  from  300  to  600  nm.  This 
measurement  was  made  in  a  synchronous  scan  mode,  whereby 
the  excitation  and  emission  gratings  are  moved  simultaneously. 
The  diffuse  reflectance  measurements  were  made  at  an  incre¬ 
ment  of  5  nm,  with  a  signal  acquisition  time  of  O.ls/wavelength 
and  at  a  PMT  high  voltage  of  340  or  450  V.  The  difference  in  the 
PMT  high  voltage  is  due  to  two  different  instruments  being  used 
in  the  course  of  the  study  (both  identical  models),  each  having 
slightly  different  throughputs.  The  instruments  were  otherwise 
identical  and  were  calibrated,  such  that  this  difference  does  not 
introduce  artifacts  into  the  data.  The  diffuse  reflectance  spec¬ 
trum  was  corrected  for  the  wavelength-dependent  system  re¬ 
sponse  and  the  throughput  of  the  instrument  by  normalizing  it  to 
that  of  a  Spectralon  99%  reflectance  puck  (SRS-99-010,  Lab- 
sphere,  Inc.,  North  Sutton,  NH),  which  was  measured  with  the 
probe  in  contact  with  the  puck  (no  coupling  media  was  used). 


Kfomial  fibt'ouy  tissue _ 

Adctiogig _ 

Rep<a'’ative  Chniiges _ 

Fibrocystic  disease _ 

Fibroad^ioiiia _ 

Cystic _ 

Total  Normal  Adipose  Tissues 


After  each  measurement,  the  exact  site  where  the  probe  was 
positiotfed  on  each  rissue  sample  was  inked  (TMD-BK,  Triangle 
Biopftical  Sci^ps,  Durham,  NC)  and  then  cut  and  stained 
for  histdpathql^)'!  Microscopic  evaluation  of  each  histological 
^section  was  ^eif ormed  (I^,  FX)  and  consensus  diagnosis 
'!;was  reached;  When  a  saniple  exhibited  a  heterogeneous  diag- 
nosis  at  the  site  of  ni^&urement,  the  mdst  se^e  diagnosis  was 
used  (e*g.,  for  samples  which  had  both  nomml  glandular  tissue 
asjli^l  as  mali^ant  tissue,  the  sample%iag|iosis  was  desig¬ 
ned  as  bein^  malignant).  In  samples  in  wlliek  both  normal 
l||dipose  and  ffiroglandular  tissuesTwferEi^pent;^^  histology 
;Fwas  desigiftited  to  be  that  of  the  Jfe^o^^fetissue  %pe  at  the 
'Site  of  ndeasurement.  Table  I  shops  the  histolbj^cal  breakdown 
oKthe  56  breast  tissue  samples  i|enined  in  fee  optical  spec¬ 
troscopy  study.  Eleven  of  tte  malipimt  sampfis  had  a  histolog- 
^al  normal  counterpaitffOfethe  saiftep^ent. 


■■"  C  Data  Analysis 

Fig.  1  shows  |lfie  fl^;ghart  oppie  algorithm  that  was  de¬ 
veloped  to  analyze  bre^ltjissudf^uorescence  and  diffuse  re- 
flect&^/JP^cSalTwo  different  types  of  input  data  sets  were 
In  efirst  the  original  spectra  were  used.  The  pri- 
^Inary  sour^S  of  variance  in  the  original  spectra  are  the  intensity 
and  spectT§4ine  shape^jn  the  second  case,  each  fluorescence 
^Spectrurfwas  normalized  to  its  peak  intensity,  or  in  the  case 
^S|i||ie  diffuse  reflectance  spectrum,  to  its  integrated  intensity. 
Tn&jnmary  source  of  variance  in  this  case  was  attributed  to  the 
specnEl  line  shape.  The  benefit  of  preprocessing  by  normaliza¬ 
tion  is  thp  it  cancels  out  inter-patient  variations  in  intensity  and 
allows  a  direct  comparison  of  spectral  line  shape. 

Principal  Component  Analysis  (PCA)  was  employed  as  a  data 
reduction  technique  [22].  PCA  characterizes  a  majority  of  the 
variance  while  greatly  reducing  the  input  data  set  into  a  few 
orthogonal  variables.  The  principal  components  (PCs)  are  ex¬ 
tracted  such  that  the  first  principal  component  (PCI)  accounts 
for  the  largest  amount  of  the  total  variance  of  the  input  data. 
The  second  PC  (PC2)  accounts  for  the  second  largest  amount  of 
the  variance  while  being  orthogonal  to  PCI,  and  so  on.  Thereby 
PCA  projects  the  set  of  observed  data  onto  a  subspace  expanded 
by  the  PCs.  There  are  two  advantages  of  this  transformation:  1) 
the  input  data  can  be  represented  by  a  few  subsets  of  PCs  with 
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Fig.  1.  Flow  chart  of  the  algorithm  that  was  developed  to  analyze  breast 
tissue  fluorescence  and  diffuse  reflectance  spectra  (PC A — principal  component 
analysis;  SVM — support  vector  machine). 


minimal  mean  square  error,  which  reduces  the  dimensionality 
of  the  data  set;  2)  the  projection  onto  the  PC  subspace  maxi¬ 
mizes  the  separation  of  data  clusters.  PCA  can  be  carried  out  on 
concatenated  spectra  (fluorescence  emission  spectra  at  all  nine 
excitation  wavelengths  and  diffuse  reflectance  spectra).  Alter¬ 
natively,  PCA  can  be  performed  on  all  fluorescence  emission 
spectra,  one  excitation  wavelength  at  a  time  and  individually” 
on  the  diffuse  reflectance  spectra.  Both  approaches  will  vi||d 
the  same  set  of  eigenvectors  for  each  spectrum.  However,  ®1ie 
first  case,  one  set  of  PC  scores  will  be  generated  for  all  s^^tra 
per  sample,  whereas  in  the  second  case,  one  set  of  PCs  will  Jbe 
generated  for  each  spectrum  per  sample.  In  this  study^  PCA  wa^ 
performed  on  all  fluorescence  emission  spectra,  btie  excitation  " 
wavelength  at  a  time  and  separately  on  the  dyfase  reflecj^ce 
spectra  to  yield  a  set  of  PCs  for  each  spectral  data  set.  Th&sub- 
sets  of  PCs  that  account  for  the  majority  (9S%.  pr  99%)  of  the 
variance,  rather  than  all  of  the  PCs,  were  retyped  Jor  further 
processing.  Using  more  PCs  (99%  versus  the  95%  Criterion)  re-'l 
tains  a  greater  proportion  of  the  variance  of  the  input  data  set 

A  Wilcoxon  rank-sum  (unpaired)  or  rank-si^ed  (pSred)^st 
[23]  was  then  employed  to  determine  which  PCskJores  showed 
statistically  significant  differences  between  malignant  and  non- 
malignant  tissue  types.  This  particular  test  was  used^ because  it 
is  a  nonparametric  test,  which  is  appropriate  for  small  sample 
sizes  for  which  a  normal  distribution  cannot  be  assumed.  Using': 
a  standard  criterion  of  p  <  0.05,  PCs  from  all  nine  fluorescence'"^ 
spectra  and  the  diffuse  reflectance  spectrum  showed  signifi¬ 
cant  differences  between  malignant  and  nonmalignant  tissues. 
However,  in  a  clinical  setting,  it  would  be  desirable  to  limit  the 
number  of  measurement  parameters  (and  thus  the  measurement 
time)  required  for  diagnosis.  Thus,  PCs  that  showed  differences 
between  malignant  and  nonmalignant  tissue  types  at  varying 
significance  levels  of  p  <  0.05,  p  <  0.005  and  p  <  0.0005 
were  retained  as  inputs  into  the  classification  scheme.  This  was 
carried  out  to  evaluate  the  classification  accuracy  as  a  function 
of  the  number  of  PCs  (or  measurement  parameters)  retained. 

A  support  vector  machine  (SVM)  algorithm  [24]  was  used 
as  the  classification  technique.  SVM  is  a  classification  algo¬ 
rithm  based  on  statistical  learning  theory.  The  principal  idea  of 
an  SVM  is  to  determine  an  optimal  separating  hyperplane  that 
maximizes  the  margin  between  two  classes  in  a  multidimen¬ 
sional  data  space.  With  the  largest  separation  of  the  two  data 


clusters,  the  SVM  classifier  gives  a  lower  expected  risk,  which 
means  that  future  error  can  be  minimized  if  more  data  is  added 
to  the  sample  pool.  Specifically,  a  linear  kernel  SVM  classifier 
was  employed  in  this  study,  and  the  classification  was  carried 
out  on  PCs  retained  from  both  fluorescence  spectra  and  the  dif¬ 
fuse  reflectance  spectrum  or  the  diffuse  reflectance  spectrum 
only. 

A  full  cross-validation  of  the  algorithm  was  performed  in 
the  following  manner.  First,  one  sample  was  removed  from  the 
training  dataset.  The  spectra  of  the  remaining  samples  under¬ 
went  the  steps  of  preprocessing,  PCA,  Wilcoxon  testing  and 
SVM  classification.  The  optimal  PC  scores  (and  corresponding 
eigenvectors)  were  identified  from  the  Wilcoxon  test  and  the  op¬ 
timal  decision  hyperplane  was  determined  from  the  SVM  anal¬ 
ysis.  Next,  the  optimized  algorithm  was  applied  to  the  sample 
that  had  been  left  out.  In  particular,  the  PC  scores  were  calcu¬ 
lated  from  the  sample  spectra  using  the  optimal  eigenvectors 
retained  froiftthe  training  set.  Then,  the  sample  was  classified 
as  malignSnt  or  nonmalignant  based  on  which  side  of  the  op¬ 
timal  decision  hyperplane  its  PC  scores  were  distributed.  This 
method- allowed  Tor  an  unbiased  evaluation  of  the  algorithm’s 
performance  wjhen  the  sample  size  is  too  small  to  separate  into 
jndependent^ftaining  and  testing  data  sets^  This  procedure  was 
carried  all  of  th^;f6  samples.  Por  each  testing  sample, 

the  algorithm  was  trairieff  on  the  remaining  55  training  samples. 

In  the  original  sta^  of  algorithm  deveii^^ent,  a  two-step  al- 
gorifhhi  was  employed  to  differentiate  between  malignant  and 

Jmmalignant  (nbfmal  and  benign)  tissues.  involved  first 
scriminatin§  between  fibrous/gl^dia^rand  g^ippse  tissues 
and  then  d||criminating  betweei^^aligri^tiand  fibrous/glan- 
.^ular  tissues.  This  step-wise  app&ch  resulte^pin  samples  being 
Mscl^sified  at  each  step.  It  was^bserved  thala  greater  number 
oPsOTples  were  misclassified  uSg  this  ap^bach  relative  to 
that  in  which  the  algorithm-was  uSed;ta;dfectly  discriminate 
;"fetween  the  two  class^f  intefi^&e.,  ril^i;^ant  and  nonmalig- 
bnant  (normal  and  beitign)  tissue. "ipferefore,  the  latter  approach 
was  employed  inAi&tudy.  | 

It  is  desirable  lb  feldtpthe  PCs  tliitt  contribute  to  optimal  clas¬ 
sification  of  4ltfignant  nonMalignant  breast  tissues  to  fea- 
turefSAe  COTegponding  spectrum.  This  can  be  achieved  by 
fbr^rojeb&g  thel^s  back  into  the  spectral  data  space.  Re-pro- 
^&tion  of  mb  PC  of  %given  sample  into  the  spectral  data  space 
can  be  achieved  by  ffi^Jftiplying  the  PC  score  of  that  sample 
::/by  the  ^biresponding  eigenvector.  A  linear  combination  of  the 
Terpro^ted  PCs  that  account  for  100%  of  the  variance  can  faith¬ 
fully  reconstruct  the  spectrum  of  that  sample.  Since  the  objec¬ 
tive  ls:fedetenrdne  which  spectral  features  are  useful  in  the  clas¬ 
sification;  only  those  PCs  that  showed  the  statistically  most  sig¬ 
nificant  differences  between  malignant  and  nonmalignant  tis¬ 
sues  were  used  in  the  re-projection. 

in.  Results 

Fig.  2  shows  the  average  fluorescence  excitation-emission 
matrices  (EEMs)  for  malignant  (20)  [Fig.  2(a)],  normal/benign 
fibrous  (15)  [Fig.  2(b)],  and  adipose  tissues  (21)  [Fig.  2(c)]. 
Each  EEM  is  shown  on  a  log  contour  scale,  where  each  con¬ 
tour  corresponds  to  levels  of  equal  fluorescence  intensity.  The 
useful  fluorescence  information  is  to  the  right  of  the  Rayleigh 
scattering  line  seen  in  the  upper  left  quadrant  of  the  plot.  This 
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Fig.  2.  Average  fluorescence  EEMs  for  (a)  malignant  (n  =  20),  (b) 
normal/benign  fibrous  (n  =  15),  and  (c)  adipose  tissues  {n  =  21).  Each  EEM 
is  shown  on  a  log  contour  scale,  where  each  contour  corresponds  to  levels  of 
equal  fluorescence  intensity.  The  useful  fluorescence  information  is  to  the  right 
of  the  Rayleigh  scattering  line  seen  in  the  upper  left  quadrant  of  the  plot. 

plot  allows  for  the  display  of  fluorescence  intensities  acquired 
at  multiple  excitation  wavelengths  where,  the  fluorescence  spec¬ 
trum  at  each  excitation  wavelength  is  corresponds  to  a  horizon- 


Fig.  3.^  Average  diAhMYeflectance  spectra  with  standard  deviations,  for  the 
maliflant,  normal/b^rgn  fibrous  and  adipose  tissues. 

J:-  ^^4. 

^•  tal  slice  thnjwgb  the  plfot.  Four  peaks  visible  in  the  malig- 

nant  [IQg.  fibrdwsftissues  [Fig.  2(b)], 

whicli  Bccur  in  simlar  locations,  i.e.,  rflixcitation,  emission 

S^^ength  pmrs  ^(300, 340),  (340, 390)1  and  (440, 

)  nm.  Thq'ajdipose  tissue  EEM  pFig;-g|^).]  on  tife  pther  hand, 
iJiows  disrfpf  differences,  relativ#to  tha^pervediii  Fig.  2(a) 
mid  (b).  jf'ptably,  the  peak  at  (310;,  390)  nnl^is:  weakly  present 
ffiH  thd^ak  at  (360, 460)  nm  halfeen  shifted  to  approximately 
(360;  SiO)  nm.  : 

Fig.  3  shows  the  average  ditese  refleeHg#  spectra  with  stan- 
fd  deviations  for  thfemalighmiti  normal/benign  fibrous  and 
|f,  adipose  tissues.  It  cafi  jbe  seen  thftthe  average  malignant  and 
normal/benign  fi^usrti^sue  spectr|jiave  similar  intensities,  but 
subtle  differenc^nn  sp^tral  line  ilape  in  the  wavelength  re¬ 
gions  betweert3j&-350,  400-45||and  525-575  nm.  The  average 
adj^lje  tissue  spectrum  shows  decreased  intensity  at  all  wave- 
Jjenfths,  and>a  disti|ictly  different  line  shape  between  425-525 
pm  compafod  to  the  average  malignant  and  normal/benign  fi- 
r  )  brous  tiss^^rspectra.  Sote  also  that  the  variability  in  adipose 
tissue  spectra  is  significantly  less  than  that  of  fibrous  and  ma¬ 
tt  tissues. 

ifele  II  shows  the  results  of  PCA,  followed  by  a  Wilcoxon 
rank-^ip,  test  for  original  [Table  n(a)]  and  normalized 
[Table  fluorescence  and  diffuse  reflectance  spectra. 

The  PCs  flat  display  differences  between  the  malignant  and 
nonmalignant  tissues  at  different  significance  levels  are  listed. 
Each  table  is  arranged  to  show  which  PCs  meet  a  particular 
criterion  but  do  not  meet  a  higher  criterion  [for  example  in 
Table  11(b),  PC2  at  320-nm  excitation  meets  both  the  p  <  0.05 
and  p  <  0.005  criterion  but  not  p  <  0.0005  criterion,  and  so 
is  shown  in  the  p  <  0.005  column].  The  PCs  are  organized 
based  on  which  spectra  they  were  derived  from.  For  a  particular 
spectrum,  the  PC  number  and  the  variance  described  by  that 
PC  (in  parentheses)  are  indicated.  Note  that  Table  n  shows 
the  results  of  PCA  and  the  Wilcoxon  rank-sum  test  for  all  56 
samples. 
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TABLE  II 

The  Results  of  PCA,  Followed  by  a  Two-Tailed  Wilcoxon  Rank-Sum 
Test  for  (a)  Original  and  (b)  Normalized  Fluorescence  and  Diffuse 
Reflectance  Spectra.  The  PCs  That  Display  Differences  Between  the 
Malignant  and  Nonmalignant  Tissues  at  Different  Signihcance 
Levels  Are  Listed.  Each  Table  Is  Arranged  to  Show  Which  PCs  Meet 
A  Particular  Criterion  but  do  not  Meet  a  Higher  Criterion  (for 
Example  in  (b),  PC2  at  320  nm  Excitation  Meets  Both  the  p  <  0.05 
AND  p  <  0.005  Criterion  but  not  p  <  0.0005  Criterion,  and  so  Is 
Shown  in  the  p  <  0.005  Column).  The  PCs  Are  Organized  Based  on 
Which  Spectra  They  Were  Derived  From.  For  a  Particular  Spectrum, 
THE  PC  Number  and  the  Variance  Described  by  That  PC  (in 
Parentheses)  Are  Indicated.  Note  That  Table  II  Shows  the  Results  of 
PCA  AND  THE  Wilcoxon  Rank-Sum  Test  for  all  56  Samples 


3(K)imi£xc.  PC2(Q.C^>) _ PCI  (0-93) 

320muBxc.  PCI  (0.90) _ 

340muBxc.  PCI  (0.90) _ 


360imiBxc.  PCI  ((0.S9) 


cant  than  the  PC  retained  from  the  diffuse  reflectance  spectrum. 
In  Table  11(b),  there  are  a  total  of  seven  PCs  retained  from  flu¬ 
orescence  spectra  at  several  excitation  wavelengths  that  are  sta¬ 
tistically  more  significant  than  the  two  PCs  retained  from  the 
diffuse  reflectance  spectrum. 

In  addition  to  the  unpaired  Wilcoxon  rank-sum  text,  a 
Wilcoxon  rank-signed  test  of  normal/malignant  PC  pairs  within 
the  same  patient  was  performed  for  a  total  of  1 1  normal/ma¬ 
lignant  tissue  pairs.  The  Wilcoxon  rank-signed  test  resulted  in 
similar  PCs  to  that  generated  from  the  Wilcoxon  rank-sum  test. 
Specifically,  11  of  the  15  PCs  identified  as  being  significant 
from  an  unpaired  analysis  of  normalized  spectra  [Table  11(b)] 
were  also  identified  as  being  significant  from  a  paired  analysis 
of  the  normalized  spectra.  The  consistency  between  the  paired 
and  unpaired  analysis  confirms  that  the  observed  spectral 
differences  between  malignant  and  nonmalignant  tissues  are 


380  mu  £xc. 
400imiBxc. 
A20im.Bxc. 
440mnExc. 
460  nm  Bxc. 
Reflectance 


Spectra 
300  im\  Hxc. 
320  mn  Exc. 
340  mu  Exc 
360  mn  Exc. 
380  im\  Exc. 
4CK)  mil  Exc, 
420mxtExc 
440  mu  Exc. 

460  mu  Exc. 
Reflectance 


PCI  (0.9?) 
PCI  (0.991) 
PCI  (0.995) 
PCI  (0.996) 
PCI  (0.997) 


PCI  (0.72) 
PCI  (095) 
PCI  (0.75) 
PCI  (PC?) 

PC3ffl03) 


PCS  (0.12), 
PC2  (0.37) 


PCI  (0.80) 


PC2(().n) 


PCI  ^.68) 

"  PC4  "(bs>r 

PC!  (0.67) 


p<00005 
PC2  (0.12) 


PC4  (0.02) 
PCI  (0.66) 


PCS  (0,05) 


A  comparison  of  Table  11(a)  and  (b)  indi^tes  thaTfar  nidrb 
PCs  at  higher  significance  levels  (p  <  O.OOSHapdp  <  6^0005) 


are  generated  using  normalized  versus  the  origihi^|luorescence 
spectra.  Thus,  this  preprocessing  method  was  utilif^in  all  sub¬ 
sequent  analysis  of  fluorescence  spectra.  It  is  interesting  to  note 
that  both  the  original  and  normalized  fluorescence  spectra  at 
300-nm  excitation  yield  PCs  that  meet  the  p  <  0.0005  crite-^ 
rion.  In  Table  n(b),  it  can  be  seen  that  as  the  significance  level 
increases,  the  PCs  that  are  retained  from  the  fluorescence  spec¬ 
tral  data  account  for  a  smaller  percent  of  the  total  variance.  At 
a  significance  level  of  p  <  0.005,  the  most  useful  fluorescence 
spectra  are  those  measured  at  300,  320, 400, 420, 440,  and  460 
nm  excitation.  In  the  case  of  the  diffuse  reflectance  measure¬ 
ments,  the  original  spectrum  [Table  n(a)]  yields  only  one  PC 
that  meets  thep  <  0,005  criterion,  while  the  normalized  spec¬ 
trum  [Table  n(b)]  yields  two  PCs  that  meet  the  p  <  0.05  crite¬ 
rion,  but  not  the  p  <  0.005  criterion.  Thus,  both  preprocessing 
methods  were  thus  utilized  in  subsequent  analysis  of  diffuse  re¬ 
flectance  spectra. 

PCs  at  higher  significance  levels  are  generated  from  fluo- 


similar  within  and  between  patients. 

Table  in^tiows  the  classification  results  obtained  from  SVM 
analysis  qI  fee  diagnostically  important  PCs  derived  from  nor- 
malized%ectra.  The  PCs  used  are  similar  to  those  indicated  in 
TablCi  llfb).  However,  it  should  be  noted  that  Table  11  shows  the 
resiS^|)f  PCi^^pwl  the  Wilcoxon  rank-sum  test  for  all  56  sam¬ 
ples,  Wpile  fe  ffie  actual  al|orithm  development,  only  55  sam- 
i  pies  weire  used  at  a  time  in  the  trainingfi^t  (one  testing  sample 
-  was  remove^  sequentiffly  for  the  purples  of  cross-validation). 
Thus,  the  PC!s  show  in  Table  11  could  vi^^lightly  from  those 
idejt0ed  from  fe^  56  independent  traininf  ^||ls.  The  classifica- 
tidh  fate,  sensitivity  and  specificity  for  the  tftttog  and  testing 
Mzu  sets  are  shown  for  PCs  retained aff||gressryeiy  increasing 
levels  of  silhificance  ranging  fr^  p  to^:  <  0.0005. 

V  As  seen  -in  Table  n(b),  the  PCs||etained  anfe^  p  <  0.05  sig- 
nifilcaride  level  originate  from  bdifc|the  fluorel^nce  and  diffuse 
reflect^ce  spectra,  while  those  feifemed  at  hipier  significance 

Jels  originate  only  frpm^fee'flpofei^enq^  Ip^  The  PCs  in 

first  three  rows  wdi^  identified  frdnfS  set  of  PCs  that  de¬ 
scribe  95%  of  the  vai^ce  of  the^put  data  set.  The  PCs  in  the 
fourth  row  were  iden^pd  from  a  fet  of  PCs  that  describe  99% 
of  the  variance  df  the  iRpat  data  set  increasing  the  percent  vari¬ 
ance  from  95%  lb  99%  in  two  to  three  more  PCs  at  the 
p  ^  Q.0005  level*  The  training  data  set  shows  a  classification  ac- 

J^Scy  (cl^sificMpn  rate,  sensitivity  and  specificity)  that  cor- 
sponds  ^  the  avefa|e,  and  standard  deviation  of  that  obtained 
from  trairjpg  with  a  of  56  data  sets,  each  with  a  total  of  55 
^|mple^|bne  sample  at  a  time  of  a  total  of  56  samples  was  left 
^||fqf  &  purposes  of  cross-validation).  The  classification  ac- 
ciMg|,  of  the  testing  data  set  reflects  the  unbiased  performance 
of  tK%.|lgorithm  on  56  samples  through  cross-validation. 

Tabfa|n  shows  that  cross-validation  of  the  algorithm  results 
in  a  5%-20%  decrease  in  the  classification  accuracy.  It  is  inter¬ 
esting  to  note  that  that  the  classification  accuracy  of  the  testing 
data  set  shows  a  smaller  deviation  relative  to  that  of  the  training 
data  set  as  the  number  of  PCs  used  in  the  SVM  analysis  is  re¬ 
duced  (by  increasing  the  level  of  significance).  Furthermore,  the 
classification  accuracy  of  testing  data  set  increases  slightly  for 
increasing  levels  of  significance.  Selecting  the  diagnostically 
relevant  PCs  at  the  p  <  0.0005  level  of  significance  from  a  set 
of  PCs  that  describe  99%  as  opposed  to  95%  of  the  variance  of 


rescence  spectra  than  from  the  diffuse  reflectance  spectra.  In  the  input  data  set  gives  rise  to  an  increase  in  the  number  of  PCs 
Table  n(a),  there  is  one  PC  retained  from  the  fluorescence  spec-  used  for  SVM  classification.  This  results  in  an  improvement  of 
trum  at  300  nm,  excitation,  which  is  statistically  more  signifi-  the  overall  classification  accuracy,  as  well  as  the  smallest  dif- 
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TABLE  III 

The  Classification  Results  Obtained  From  SVM  Analysis  of  the  Diagnostically  Important  PCs  Derived  From  Normalized  Spectra.  The 
Classihcation  Rate,  Sensitivity  and  Specificity  for  the  Training  and  Testing  Data  Sets  Are  Shown  for  PCs  Retained  at  Progressively 
Increasing  Levels  of  Significance  Ranging  From  p  <  0.05  to  7?  <  0.0005.  As  Seen  in  Table  11(b), the  PCs  Retained  at  the  p  <  0.05  Signihcance 
Level  Originate  From  Both  the  Fluorescence  and  Diffuse  Reflectance  Spectra,  While  Those  Retained  at  Higher  Signihcance  Levels 
Originate  Only  From  the  Fluorescence  Spectra.  The  PCs  in  the  First  Three  Rows  Were  Identified  From  a  set  of  PCs  That  Describe  95%  of  the 
Variance  of  the  Input  Data  Set.  The  PCs  in  the  Fourth  Row  (Asterisks)  Were  Identihed  From  a  set  of  PCs  That  Describe  99%  of  the  Variance 
OF  THE  Input  Data  Set.  The  Training  Data  Set  Shows  a  Classihcation  Accuracy  (Classification  Rate,  Sensitivity  and  Specificity)  That 
Corresponds  to  the  Average  and  Standard  Deviation  of  That  Obtained  From  Training  With  a  Total  of  56  Data  Sets,  Each  With  a  Total  of  55 
Samples  (One  Sample  at  a  Time  of  a  Total  of  56  Samples  was  Left  out  for  the  Purposes  of  Cross-Validation).  The  Classification  Accuracy  of 
THE  Testing  Data  Set  Reflects  the  Unbiased  Performance  of  the  Algorithm  on  56  Samples  Through  Cross-Validation 


Sigtiificcoice  Level 

Number  of  PCs 

Ckissifkatian  Rate  (®/o) 

Sensitivity  (%) 

Specificity  (%) 

Training  Data 

Testing  Data 

Training  Data 

Testing  Data 

Training  Data 

Testing  Data 

P<0.05 

14-15 

91.53+276 

75 

84.16±5.66 

65 

95.60±210 

80.56 

P<Q.005 

7-8 

86.95±1.87 

76.79 

77.36±4.65 

65 

92.27±3.06 

83.33 

P<0.0005 

3-4 

83.51+1.89 

79 

76.07±5.61 

70 

87.63±2.63 

83.33 

*P<0.0005 

6-7 

87.56±213 

83.93 

76.35±5.54 

70 

93.78±1.92 

91.67 

TABLE  IV 

Classircation  Results  Obtained  From  SVM  Analysis  of  the  Diagnostically  Important/PCs  Derived  From  Both  the  Original  and  Normalized 

Diffuse  Reflectance  Spectra.  This  a  Similar  ^orS^at  to  Table  III 


Pi^-procesKii^ 

Signiliccmce  Level 

Number  of 

Clas.s]fieati0n  Rate  (%) 

Sensitivity 

Specificity  ^/«) 

PCs 

Training  Data 

Testing  Data 

TrainingData 

Testing  Data 

TrainingData 

Testmg  DaAa 

Original 

Normalized 

P<0.i 

P<0.05 

2 

2 

66.62+251 

69.46±5.72 

60.71 

50 

45.43±8.85 

5L38+16.3 

30 

20 

78.3S±1.73 

79.41+5.07 

77.78 

66.67 

ference  in  the  classification  accuracy  (in  particular,  the  speci^^r  ized  fluorescence  spej(^a'at  300  nm,  efe^tijdn),  for  malignant 
ficity)  between  the  training  and  testing  data  sets.  A  comparison  and  ndnmali^^^t  tipties,  and  the  difFef6hc|’spectra  (between 
of  the  classification  accuracy  of  PCs  retained  at  the  p  <  QfPOOS  maUgnant  and  nl**^ — - *  An^\i 


;  tissues)  forSacE^ase  [Fig.  4(b)]. 


and  p  <  0.05  level  of  significance  in  Table  HI  indicates  that  ttgspectral  representations  are  re-projections  from  the  PC  sub¬ 
combining  the  normalized  fluorescence  and  diffuse  reflectA^e  space  onto  normalized  spectral.data  space.  They  are  calcu- 
spectra  does  not  improve  the  classification  accuraegof  an  alg^^i;  dated  by  skifiiy  multiplying  the  PJi"  scores  fhj  thil  case,  PC2) 
rithm  based  on  normalized  fluorescence  spectra  afcfie.  Mor  all  spnples  by  the  associatS^  eigenvectoJ;  and  then  aver- 

Table  IV  shows  classification  results  obtai&ed  from  SjVM  aging  spectral  representations  for  each  tipne  type  (malig- 
analysis  of  the  diagnostically  important  PCsfeived  frorfSbth  naift^apd  nonmalignant).  The  difference  spe^lnn  was  obtained 
the  original  and  normalized  diffuse  reflecfhnge,  spectral  This  tjy  suWacting  (point-b^^0ih^:^e  avefageWhmalignant  spec- 
table  has  a  similar  format  to  Table  in.  Again ^ed^fipselected  item  from  the  averagetfcalignkfif^iectrum.  PC2  was  selected 
for  this  analysis  are  similar  to  those  slli^n  in  Tafeleit.  Table  for  the  spectral  reprelintation  beS^e  the  Wilcoxon  rank-sum 
shows  that  only  one  PC  was  retained  frbTn:tbe  original  diffuse  test  indicated  thafPGZ  displays  tfe  statistically  most  signifi- 
reflectance  spectra  [Table  11(a)]  and  two^BEs  were'ml^ed  cant  differences^fetween  maUgn^t  and  nonmalignant  tissues. 

It  should  be  pbinted  out  thaf  ^though  PCI  accounts  for  more 
of  |h|  vaid^ce  fiflhe  spectral  data  than  PC2,  it  also  reflects  the 
yBance  p|||ent  m  both  malignant  and  nonmalignant  tissues, 
i  In  Fig.  4|ft|,  the  aVnmge  normalized  fluorescence  spectrum  of 
diagnostically  relevant  PCs  from  the  original  diffuse  reflectance^  the  malignant  tissues  shbws  a  red  shift  in  the  fluorescence  peak 


from  the  normalized  diffuse  reflectance  speb^';[Table"^ 
when  a  criterion  of  p  <  0.05  was  used.  In  orefer  to  carry  out 
SVM  analysis,  at  least  two  features  (PCs)  are  ne^i^ed.  Thus, 
the  p  <  0.1  criterion  was  used  in  order  to  select  af  least  two 


spectra.  ^ 

Table  IV  shows  that  cross-validation  of  the  algorithm  results 
in  a  10%-30%  decrease  in  the  classification  accuracy.  The  clas¬ 
sification  accuracy  of  the  testing  data  set  shows  a  smaller  devia¬ 
tion  relative  to  that  of  the  training  data  set  for  PCs  derived  from 
the  original  diffuse  reflectance  spectra.  In  addition,  the  classi¬ 
fication  accuracy  of  testing  data  set  for  PCs  extracted  from  the 
original  spectra  is  higher  than  that  for  PCs  derived  from  the  nor¬ 
malized  spectra.  A  comparison  of  the  classification  accuracy  of 
PCs  retained  at  the  p  <  0.05  or  p  <  0.1  level  of  significance  in 
Table  IV  and  p  <  0.0005  level  of  significance  in  Table  El  indi¬ 
cates  that  using  normalized  fluorescence  spectra  yields  signifi¬ 
cantly  higher  classification  accuracy  than  using  either  the  orig¬ 
inal  or  normalized  diffuse  reflectance  spectra  alone. 

Fig.  4  shows  the  average  normalized  fluorescence  spectra  at 
300-nm  excitation  [Fig.  4(a)]  and  the  average  spectral  represen¬ 
tations  reflected  by  PC2  (which  was  obtained  from  the  normal- 


^d  a  lower  intensity  between  450-550  nm  relative  to  nonma¬ 
lignant  tissues,  These  differences  are  reflected  in  the  average 
spectel  representation  of  the  malignant  and  nonmalignant  tis¬ 
sues.  Thg, differences  seen  in  the  average  malignant  and  nonma¬ 
lignant  Slue  spectral  representations  track  the  differences  ob¬ 
served  in  their  corresponding  average  normalized  fluorescence 
spectra,  but  do  not  necessarily  replicate  the  actual  spectra  it¬ 
self,  since  only  one  PC  is  used  in  the  re-projection.  Fig.  4(b) 
shows  the  difference  spectra  (between  malignant  and  nonma¬ 
lignant  tissues)  for  both  the  average  normalized  spectra  and  av¬ 
erage  spectral  representations  shown  in  Fig.  4(a).  It  can  be  seen 
that  the  two  difference  spectra  are  very  similar  to  each  other, 
as  expected.  This  verifies  that  the  PC2  displays  the  most  sig¬ 
nificant  difference  between  the  two  tissue  types.  The  difference 
spectra  indicate  that  the  most  significant  differences  between 
malignant  and  nonmalignant  tissues  are  observed  near  emission 
wavelengths  of  350  and  500  nm,  when  an  excitation  wavelength 
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•  - '  -  Average  spectral  representation  of  PC2  for  malignant  samples 

—  Average  normalized  spectra  of  mafignant  samples 

—  Average  spectral  representation  of  PC2  for  non-malignant  samples 

Average  normalized  spectra  of  non-malignant  samples _ 


yy' . '-O- 


400  450  500 

Emission  Wavelength  (nm) 


Differertce  spectra  obteined  from  normalized  spectra 
Difference  spectra  obtaJned  from  spectral  representation  of  PC2 


350  400  450  500  650  ^ 

Emission  Wavelength  (nm)  I 

Fig.  4.  (a)  The  average  normalized  fluorescence  spe^^jt^OO  nm^i^xciMtion, 
and  the  average  spectral  representations  reflected  by  PC2^|dch  was  bifened 
from  the  normalized  fluorescence  spectra  at  300  nm,  excitSi^j|[i^  for  malignant 
and  nonmalignant  tissues,  and  (b)  the  difference  spectra  (b^^^n  malignant 
and  nonmalignant  tissues)  for  each  case.  The  spectral  rep^&itations  are 
re-projections  from  the  PC  subspace  onto  the  normalized  spectral  data  space. 
They  are  calculated  by  simply  multiplying  the  PC  scores  (in  this  case,  PC2j 
for  all  samples  by  the  associated  eigenvector  and  then  averaging  the  spectr3f^ 
representations  for  each  tissue  type  (malignant  and  nonmalignant).  The 
difference  spectrum  was  obtained  by  subtracting  (point-by-point)  the  average 
nonmalignant  spectrum  from  the  average  malignant  spectrum. 


of 300  nm  is  used.  Spectral  representations  obtained  from  PC4  at 
400-nm  excitation,  and  PCS  at  460-nm  excitation  are  not  shown 
here,  but  indicate  that  the  maximal  differences  between  malig¬ 
nant  and  nonmalignant  tissues  are  observed  at  emission  wave¬ 
lengths  of  460  and  540  nm,  respectively. 


those  spectra  that  are  most  useful  in  discriminating  malignant 
and  nonmalignant  tissues  were  identified.  The  analysis  indi¬ 
cated  that  only  four  of  the  ten  measured  spectra  are  required  to 
maximize  classification  accuracy.  These  include  fluorescence 
spectra  at  excitation  wavelengths  of  300,  400,  420,  and  460 
nm.  This  has  important  implications  in  clinical  applications, 
where  both  speed  and  cost  are  important,  in  that  the  fewer 
wavelengths  that  have  to  be  measured,  the  faster  and  less 
complex  the  instrument  can  be. 

One  problem  with  the  particular  approach  employed  here 
for  analysis  is  that  it  can  be  difficult  to  determine  the  biologic 
basis  of  the  differences  observed  in  the  fluorescence  spectra 
of  malignant  and  nonmalignant  tissues.  However,  an  under¬ 
standing  of  the  fluorophores  that  contribute  to  the  fluorescence 
spectra  at  the  four  optimal  excitation  wavelengths  can  yield 
some  insight  into  the  sources  of  endogenous  contrast  present. 
The  excitation  wavelengths  used,  ranging  from  300  to  460  nm, 
allow  for  thfil characterization  of  a  number  of  fluorophores, 
including'^liyptophan,  NAD(P)H,  flavoproteins  and  collagen 
[17].  Bf'S^ing  a  r|gge  of  excitation  wavelengths,  one  is  able 
to  chSnacterize  jfltey  or  all  of  the  biological  fluorophores 
preStin  brealt^tissue.  The  fluorescence  excitation-emission 
wavelepths  identified  as  Jedng  diagnostic,  from  the  PCA-SVM 
^l^lgorithl^^Uggest  that  tfi^^importantiftuorophores  for  breast 
X  cancer  dia^osis  are  ^|)tophan,  colla^e%;N^^H  and  flavopro¬ 
teins.  This  quajjtatiye  analysis  does  not  fif^|dde  the  possibility 
that  the  differeriBgs  in  the  fluorescence  spedtS^f  malignant  and 
mhutialignant  tissSes  may  also  be  influenced %^e  absorption 
Z'&’f  hemoglobin. 

As  disculsed  previously,  the  fjfcescehee  spectra  ^  300  nm, 
"^^citatiofifappear  to  be  most  dia|hostic  of  breast  cancer.  Fig.  4 
sli0ws  a  red  shift  in  the  fluorescence  peak  at  3|0  nm  (tentatively 
as^^d  to  tryptophan  fluoresceriqa)^  and  a  defease  in  fluores¬ 
cence  at  the  smaller  pegSf;cK3Cunin|''n??ui®^  emission 
/  tentatively  assigned  liiS^ADH4|feorescence)  in  the  malignant 
pi'samples,  relative  to  ifenmalignartelbsues.  This  finding  is  con- 
F  sistent  with  previous  %rk  by  Yang1%  a/.  [10],  who  found  that  at 
an  excitation  wayelenglheof  300  ml  the  fluorescence  emission 
peak  at  340  rto  '^as  red-5hflK^^of  malignant  samples  relative 
normal  adipose  tissues  and  the  fluorescence  at  around 
^^^@-nm  ‘erassiori'^^_  decreased  in  malignant  relative  to  fibrous 
,1  tissues.  It  OTy  be  tlii||he  shifted  peak  is  useful  in  discriminating 
the  malighSiit  tissues  fife  adipose,  and  that  the  decreased  fluo- 
Xrescenc^  at  450  nm  is  useful  in  discriminating  malignant  tissue 
from  hfegn  fibrous  tissues,  and  so  the  combination  of  the  two 
features  makes  this  particular  PC  useful  in  discriminating  ma¬ 
lignant  and  nonmalignant  tissues. 

Diffu^  feectance  can  characterize  some  of  the  same 
optical  features  exploited  in  fluorescence  spectroscopy,  notably 
absorption  due  to  various  biological  chromophores,  and  tissue 
scattering  properties.  The  primary  absorbers  are  oxygenated 
and  deoxygenated  hemoglobin,  which  are  present  in  blood. 
These  compounds,  which  have  slightly  different  absorption 


rv.  Discussion 

Optical  spectroscopy  in  the  UV-VIS  spectral  range  is 
demonstrated  to  be  useful  in  discriminating  malignant  and 
nonmalignant  human  breast  tissues.  By  measuring  fluorescence 
spectra  at  multiple  excitation  wavelengths  as  well  as  diffuse 
reflectance  spectra  and  performing  a  separate  analysis  on  each. 


spectra,  can  be  used  as  an  indicator  of  tissue  oxygenation 
[17].  /3-carotene  is  another  chromophore  and  is  present  in 
adipose  tissues,  thus  providing  a  unique  indicator  of  this  tissue 
type.  The  primary  advantage  of  this  technique  is  that  diffuse 
reflectance  signals  are  several  orders  of  magnitude  greater 
than  the  weak  endogenous  fluorescence  signals  measured  from 
tissue.  Therefore,  if  diffuse  reflectance  measurements  were 
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capable  of  discriminating  malignant  and  nonmalignant  tissues 
at  similar  accuracies  to  that  of  fluorescence  measurements, 
there  would  be  a  considerable  increase  in  cost-effectiveness. 
However,  using  a  similar  classification  technique,  diffuse 
reflectance  spectroscopy  was  not  capable  of  discriminating 
between  malignant  and  nonmalignant  tissues  with  the  same 
accuracy  as  fluorescence  spectroscopy.  Therefore,  there  is  an 
apparent  advantage  to  using  fluorescence  in  addition  to,  or 
instead  of  diffuse  reflectance.  This  advantage  is  likely  due  to  the 
additional  chemical  specificity  that  arises  from  the  fluorescence 
characterization  of  the  many  biologic  fluorophores  intrinsically 
present  in  tissue. 

Bigio  et  al  reported  better  classification  accuracies  for  dis¬ 
criminating  between  malignant  and  nonmalignant  breast  tissues 
using  diffuse  reflectance  spectroscopy  in  a  pilot  study  [16],  as 
well  as  in  a  more  recently  published  abstract  [25].  However, 
it  is  difficult  to  compare  classification  accuracies  achieved  in 
our  study  to  that  carried  out  by  Bigio  and  co-workers  due  to 
differences  in  the  clinical  protocol  (the  studies  by  Bigio  and 
co-workers  were  carried  out  on  tissues  in  vivo),  probe  geom¬ 
etry,  spectral  range  (the  spectral  measurements  were  made  from 
330-750  nm)  and  methods  of  analysis  employed  in  the  two 
studies.  What  can  be  concluded  from  the  work  presented  here  is 
that  fluorescence  spectroscopy  provides  superior  classification 
accuracy  compared  to  diffuse  reflectance  spectroscopy  when 
evaluated  on  the  same  samples  using  the  same  methods  of  anfr 
ysis. 

A  novel  nonparametric  algorithm  was  developed  to  identify 
and  incorporate  the  optimal  spectral  features  for  discriminating 
between  malignant  and  nonmalignant  breast  tissues;  The  use  of 
the  Wilcoxon  test  and  SVM  classification  schein4  does  not  re¬ 
quire  a  priori  knowledge  about  the  sample  distriSdtion  andj|qs, 
can  be  utilized  for  the  analysis  of  small  samplS;$izes.  Crps^al- 
idation  of  the  algorithm  highlights  the  importpce  of  o^flining 
an  unbiased  estimate  of  the  algorithm’s  classiffdpqn accuracy. 
Cross-validation  techniques  may  indicate  that  tlif  n’aining  data 
set  is  not  sufficient  to  fully  characterize  tfie  variance  of  each  sub¬ 
group  (malignant  and  nonmalignant  tissues).  Thus,  by  femovmg 
one  data  point,  the  decision  hyperplane  is  altered  enough  to 
misclassify  some  additional  samples.  Indeed,  iri%xap^^  the 
misclassified  samples,  it  is  clear  that  some  unddlr^xesented 
sub-groups  of  tissue  types  are  most  frequently  mistlassified. 
There  were  five  malignant  samples  that  were  consistently  mis- . 
classified  (once  in  cross  validation  and  repeatedly  in  the  training' 
set).  Of  these  five  samples  two  were  carcinoma  in  situ  (CIS) 
and  two  were  lobular  cancers.  These  samples  represent  the  total 
number  of  CIS  and  lobular  cancers  in  the  entire  sample  pool. 
Upon  cross-validation,  one  of  these  samples  would  be  removed 
prior  to  training  the  algorithm,  leaving  just  one  other  sample  of 
this  type  to  aid  in  its  classification.  Thus  the  total  variance  ac¬ 
counted  for  by  these  types  of  cancers  is  not  well  represented 
in  the  training  set,  which  results  in  the  frequent  misclassifica- 
tion  of  these  particular  cancers.  Upon  the  inclusion  of  more  data 
points  from  these  sample  types,  the  unbiased  (testing)  classifi¬ 
cation  accuracy  may  approach  the  accuracy  achieved  with  the 
training  data  set. 

Typically,  the  optical  spectral  features  of  different  tissue  types 
are  not  fully  separable,  which  gives  rise  to  imperfect  classi¬ 
fication,  The  classification  accuracy  is  also  dependent  on  the 
algorithm  used.  Given  a  constant  sample  size,  one  approach 


to  improving  the  classification  accuracy  may  be  to  incorporate 
more  features  into  the  classifier  (such  as  PCs).  As  depicted  in 
Table  IV,  increasing  the  number  of  inputs  did  not  significantly 
improve  the  classification  accuracy.  An  alternative  approach  to 
improving  the  classification  accuracy  is  to  use  nonlinear  rather 
than  the  linear  methods  of  analysis  employed  here,  i.e.,  PCA 
and  linear  SVM.  For  example,  using  a  nonlinear  kernel  in  SVM 
classification  may  improve  the  classification  accuracy,  since  the 
nonlinear  kernel  SVM  is  expected  to  be  superior  to  the  linear 
kernel  in  dealing  with  nonseparable  cases.  This  is  beyond  the 
scope  of  the  current  study,  but  will  be  explored  in  future  inves¬ 
tigations. 

One  of  the  goals  of  future  work  is  to  cany  out  clinical  studies 
to  explore  the  effectiveness  of  optical  spectroscopy  for  detecting 
breast  cancer  during  core  needle  biopsy.  In  particular,  side-firing 
fiber-optic  probes  will  be  designed  and  fabricated  for  use  in 
a  Mammatome  core  biopsy  needle.  The  tip  of  the  needle  has 
a  side-facing  aperture  (19  mm  long  and  2  mm  wide),  through 
which  pieces  of  tissue  can  be  vacuum  suctioned  when  the  needle 
is  insert^tnto  the  breast.  The  side-firing  probes  can  be  incor¬ 
porated  into  the  to  measure  fluorescence  and  diffuse  re¬ 
flectance  specn^of  breast  tissues  after  the  needle  is  positioned 
in  the%^  qficpncem  in  the  breast,  but  before  biopsy.  The  di- 
^:agnosis  based' bn  the  sp^t&l  features  be  compared  to  his- 
tology  of  the  Gprresporiding  biopsy  td^y^aluate  the  diagnostic 
potential  of  this  techmque. 

Itj^ould  also  be  desirable  to  gain  a  quadfitMye  understanding 
of  the  biological  basis  for  the  observed  speiM^  differences, 
.4^hich  are  difficult  to  determine  dire^fiyi&om  the^ue  spectra, 

;  As  such,  Ac  use  of  systems  emplby^g  cell  fines  aSilfissue  cul- 
"Itores  arpiping  explored.  A  stepwi$e  approached  this  technology 
is'w^ii&ted,  starting  with  simpl&ystems  and  adding  additional 
levlls  of  complexity.  For  example;  a  normal  ^d  malignant  cell 
monolayer  could  first  be  ifeaJgcd  td  es^bhsli  the  intrinsic  flu- 
/drescence  properties  dt  these  types.  Then  the  same  cells 
-  could  be  grown  in  afj-D  collagW  matrix,  adding  the  effects 
of  the  interaction:^ith  the  extraceUular  matrix.  Finally,  these 
cells  could  be  implanted  into  nuddinice  models  to  explore  the 
added  effects  df  the  vasculature  ind  adipose  tissue.  Such  a  step- 
wia|J^proach  could  allow  one  to  compare  simpler,  better  un- 
d^tood  models  tO:  what  has  been  characterized  in  tissue,  and 
4l|Ould  allow'^pe  to  determine  what  tissue  constituent  is 

{hiost  resppiisible  for  fhe  observed  spectral  differences.  This  is 
'  va  promigdhg  avenue  of  research,  and  applied  to  the  breast  could 
allbw^^  an  increased  understanding  of  the  underlying  mecha- 
nisna  fqr  how  optical  spectroscopy  is  able  to  diagnose  this  dis- 
ease.^'C^? 
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